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Abstract. An important consequence of the discovery of giant magnetoresis- 
tance in metallic magnetic multilayers is a broad interest in spin dependent effects 
in electronic transport through magnetic nanostructures. An example of such 
systems are tunnel junctions — single-barrier planar junctions or more complex 
ones. In this review we present and discuss recent theoretical results on electron 
and spin transport through ferromagnetic mesoscopic junctions including two or 
more barriers. Such systems are also called ferromagnetic single-electron tran- 
sistors. We start from the situation when the central part of a device has the 
form of a magnetic (or nonmagnetic) metallic nanoparticle. Transport charac- 
teristics reveal then single-electron charging effects, including the Coulomb stair- 
case, Coulomb blockade, and Coulomb oscillations. Single-electron ferromagnetic 
transistors based on semiconductor quantum dots and large molecules (especially 
carbon nanotubes) are also considered. The main emphasis is placed on the spin 
effects due to spin-dependent tunnelling through the barriers, which gives rise to 
spin accumulation and tunnel magnetoresistance. Spin effects also occur in the 
current-voltage characteristics, (differential) conductance, shot noise, and others. 
Transport characteristics in the two limiting situations of weak and strong cou- 
pling are of particular interest. In the former case we distinguish between the 
sequential tunnelling and cotunnclling regimes. In the strong coupling regime we 
concentrate on the Kondo phenomenon, which in the case of transport through 
quantum dots or molecules leads to an enhanced conductance and to a pronounced 
zero-bias Kondo peak in the differential conductance. 
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1. Introduction 

Since last few decades one can observe a common pursuit towards miniaturization of 
electronic systems. According to the empirical Moore's law, the number of transistors 
per microchip is doubled every two years. It is however obvious that the possibility 
of further miniaturization will be stopped in the near future due to the loss of chips' 
stability when the device components achieve critical dimensions in the nanometer 
range. Thus, the most challenging task facing contemporary science and technology 
is to implement structures, alternative to the silicon-based devices, whose size could 
be reduced further. 

The same tendency can also be observed in the development of hard discs, 
whose memory cells become smaller and smaller every year. This generates new 
challenges related to information reading/writing. An important step towards further 
miniaturization was the discovery of the giant magnetoresistance (GMR) effect. The 
GMR was discovered in 1988 in artificially layered metallic structures consisting of 
ferromagnetic 3d films separated by nonmagnetic metallic layers [TJ [2J. It turned 
out that electrical resistance of magnetic metallic multilayers depends on their 
magnetic state, and usually drops when magnetic configuration varies from antiparallel 
alignment to the parallel one [H [2j El 0] . The GMR effect occurs for current flowing in 
the film plane as well as perpendicularly to it. This phenomenon turned out to be very 
useful for applications in highly sensitive read heads, and allowed reading information 
from smaller memory cells using their weak magnetic field. 

The discovery of GMR initiated broad interest in spin polarized electronic 
transport in nanoscopic systems. It turned out that electron spin provides an 
additional degree of freedom, which considerably broadens the range of applications 
of mesoscopic systems in novel electronic devices. The spin-based nanoelectronics 
- called now spin electronics or shortly spintronics - is a relatively new area of 
mesoscopic physics dealing with the interplay of charge and spin degrees of freedom 
[3 O [3 [HI OH [lOl [TTJ [12] . Although a lot of theoretical and experimental works on the 
spintronic properties of mesoscopic systems have been carried out, this field is still in 
an early stage of development. 

An effect similar to the aforementioned GMR also occurs when nonmagnetic 
metallic layer in a trilayer structure is replaced with a nonmagnetic insulating barrier, 
and the current flows owing to the phenomenon of quantum-mechanical tunnelling 
through the barrier. This effect was discovered in 1975 in ferromagnetic planar 
junctions by Julliere [13] . and is of current interest due to applications in magnetic 
storage technology (Magnetic Random Access Memories) and in spin-electronics 
devices [SJ [7J HI M E3 OH CE2] • As in the case of GMR, the tunnel magnetoresistance 
(TMR) consists in a decrease (increase is also possible) in the junction resistance 
when magnetic configuration of the junction changes from the antiparallel to parallel 
one. Tunneling in complex junctions, particularly in mesoscopic ones, where charging 
effects become important, is still not fully explored. A specific kind of such systems 
are double-barrier junctions with a small central electrode (called an island). Such 
systems are known as single-electron transistors (SETs), mainly because electrons 
in a biased device flow one by one through the system and the transfer of single 
electrons can be controlled by a gate voltage. Electronic transport in such devices 
was extensively studied in the past decade, but mainly in the nonmagnetic limit 

Recent experiments on magnetic nano-structured materials revealed new 
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phenomena associated with the interplay of ferromagnetism and discrete charging 
effects. A typical example is a ferromagnetic single-electron transistor, i.e., a small 
metallic nanoparticle (semiconducting quantum dots or molecules are also possible) 
coupled by tunnel barriers to ferromagnetic electrodes f53J [51] . First ferromagnetic 
single-electron transistors were fabricated by Ono et al. [25] [26] and later by 
Bruckl et al. [27j . Transport in ferromagnetic single-electron transistors with 
nonmagnetic metallic islands - both normal and superconducting - was also measured 
[28] [29] 130] [31] . One should bear in mind, that the interplay of spin and charge effects 
was already studied long time ago in granular systems, in which magnetic nanoparticles 
were randomly dispersed in a nonmagnetic matrix |32] . Recently granular films were 
investigated again by several groups [33] [34] [35] [36] [37] [38] [39] [40] , but both the size 
of the grains and also their location were strictly controlled. 

In the case of sufficiently large metallic islands (but still in the nanometer 
range), discrete structure of the electronic states in the grain is not resolved and 
is irrelevant. To observe the discrete electronic states in transport characteristics one 
should either diminish size of the metallic nanoparticles [41] , or use semiconducting 
quantum dots based on two-dimensional electron gas 02] 03] [44] 05] . An alternative 
strategy is to use ferromagnetic semiconducting materials instead of metallic ones 
as the electrodes [46]. Magnetic impurities in the middle of the tunnel barrier of a 
ferromagnetic tunnel junction 07[ 05] 021 [SD] also can be considered as quantum dots 
with a very strong Coulomb interaction. Another group of single-electron devices are 
molecular ferromagnetic transistors [51] and, especially, ferromagnetically contacted 
carbon nanotubes [52l [53l El] Ell ESI Ell Ell ESI (6^ 

In the following we will review basic transport characteristics of ferromagnetic 
single-electron transistors in the sequential tunnelling, cotunneling, and strong 
coupling (Kondo) regimes. In particular, we will consider such properties of the device 
like charge and spin currents, tunnel magnetoresistance, spin accumulation, shot noise, 
Kondo effect, and others. In section 2 we review basic principles of single-electron 
transport in the case of a large metallic island attached to ferromagnetic leads. The 
size of the island is however small enough so that the charging energy is the dominant 
energy scale in the system. The limits of fast and slow spin relaxation in the island 
are also discussed. Then we consider shot noise and the role of discrete electronic 
structure in such devices. In section 3 we consider devices based on double metallic 
islands. In turn, in section 4 we deal with electronic transport through single-level 
quantum dots (QDs) in the sequential tunnelling and cotunneling regimes. Transport 
through multi-level QDs, including transport through two-level quantum dots and 
carbon nanotubes attached to ferromagnetic leads is discussed in section 5. Transport 
through quantum dots in the Kondo regime is briefly addressed in section 6. Final 
conclusions are in section 7. 

1.1. Basic concepts 

The systems considered in this review consist of a mesoscopic central part (island) 
coupled by tunnel barriers to external ferromagnetic leads. The central part is 
characterized by an addition energy, which corresponds to the energy needed for 
adding a single electron, and includes contributions from the electrostatic charging 
energy and the discrete single-particle level separation. If the charging energy is 
the most relevant energy scale, the systems exhibit the single-electron charging effects 
[T5l [63 , 64, 65, 66, 67, 68, 69]. An electron can tunnel to the central part only when the 
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first-order tunneling second-order tunneling 



Figure 1. A sketch illustrating (a) a single first-order (sequential tunnelling) 
process, and (b) a second-order (cotunneling) transport process. 



energy provided by the transport voltage surpasses the corresponding charging energy; 
otherwise, the current is exponentially suppressed and the system is in the Coulomb 
blockade regime [HI E01 EI]- Once the bias voltage is larger than the threshold 
voltage, the electrons can tunnel one by one through the system leading to the step- 
like current-voltage characteristics - so-called Coulomb staircase. The blockade can 
also be overcome by applying a gate voltage V g that leads to sawtooth-like variation 
of electric current with V g - so-called Coulomb oscillations [15] [63J [6H [671 E2] • 

When the electrodes are made of a ferromagnetic material, the system exhibits 
further interesting phenomena resulting from the interplay of charge and spin degrees 
of freedom [731 El]- In particular, the tunnelling current flowing through the 
system depends on the relative alignment of the magnetic moments of ferromagnetic 
electrodes, giving rise to the TMR effect 25, 26 , which is described quantitatively by 
the ratio 

TMR ^ AP ' fiF =^, (1) 
Rp Iap 

where Rp and Rap denote the total system resistance in the parallel and antiparallel 
magnetic configurations, respectively, and Jp, Iap are the corresponding currents. A 
simple theoretical model of TMR was introduced by Julliere [T3], who considered a 
single planar ferromagnetic tunnel junction and showed that TMR in such a device is 
given by TMR Ju " = 2plPr/(1 — PlPr), where the spin polarization p r of the lead r 
(r = L, R for the left/right lead) is defined as 

Pr+ — Pr- ,„s 

Pr = , (2) 

Pr+ + Pr- 

with p r ± being the spin-dependent density of states of lead r for the spin-majority 
(+) and spin-minority (— ) electrons. 



1.2. Transport regimes 

In the following considerations we will distinguish between the three different transport 
regimes: 

Sequential tunnelling - In the regime of weak coupling between the island and 
leads, and out of the Coulomb blockade, electron transport is dominated by processes 
of the first order in the coupling parameter. Electrons flow then consecutively one by 
one due to the tunnelling events, and the transition rate from an initial state \i) to a 
final state |/) can be determined from the Fermi golden-rule, 

a i ^f = ^\(i\H T \f)\ 3 6(e i -e / ), (3) 
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where is the relevant tunnelling Hamiltonian, and Si (ef) is the energy of the 
initial (final) state. This transport regime is known as the sequential tunnelling regime 
[THUS]- An example of a single first-order process is sketched in Fig.rjju 

Cotunneling - Although the sequential tunnelling in the Coulomb blockade regime 
is exponentially suppressed, the current still flows due to higher-order tunnelling 
processes involving tunnelling of a single, two or more electrons via intermediate 
virtual states [73 [76j [77] . These processes are known as cotunneling. An exemplary 
cotunneling process is illustrated in Fig. [TJd. 

The second order perturbation theory gives the cotunneling rate from an initial 
state \i) to a final state |/) [76] 



where the summation is over all virtual states |g), and s a is the energy od the state a 
(a = i, f,q). In the Coulomb blockade regime this tunnelling rate is only algebraically 
suppressed, contrary to the sequential transport processes which are then suppressed 
exponentially. Because of that even at low temperatures and in the strong Coulomb 
blockade regime the rates of cotunneling processes do not vanish. The second-order 
corrections become also important on resonance for intermediate coupling strengths. 

Strong coupling - For strong coupling of the metallic island to electrodes, 
the tunnelling processes lead to logarithmic corrections to conductance, and the 
perturbation theory fails at the degeneracy points of two consecutive charge states. In 
the case of quantum dots additionally the Kondo effect appears at low temperatures 
(below the Kondo temperature Tk, T < Tr-) leading to an enhanced conductance in 
the linear response regime [78]. 

2. Ferromagnetic single-electron transistors based on metallic 
nanoparticles 

In this section we review spin-polarized transport in a metallic ferromagnetic single- 
electron transistor (FM SET). The device consists of a metallic nanoparticle as 
the central electrode (island), which is coupled through tunnel barriers to external 
reservoirs of spin polarized electrons. A gate voltage is attached capacitively to the 
island, which allows to control position of the corresponding Fermi level. 

Electronic transport in nonmagnetic SETs was already extensively studied in the 
past two decades [ll[T3[7S[80l[8Tl[82[8^ Recently, the 

attention was also drawn to electron tunnelling in magnetic systems [251 121)1 1751 ITU 
1501 I5T1 W2\ , which was stimulated by recent progress in nanotechnology. It has been 
shown theoretically that some qualitatively new effects may arise from the interplay of 
charging effects and spin degrees of freedom. These include, for example, oscillations 
of TMR with increasing bias voltage, spin accumulation, enhancement of TMR in the 
Coulomb blockade regime, etc. [73 EH [SH [H3 E71 EBJ. The enhancement 
of TMR in the cotunneling regime and the oscillations of TMR as a function of the 
transport voltage have also been observed experimentally [3SJ [39J [551 HOP] • 

The capacitance C of few-nanometer-size particles is of the order of 10 -18 F 
[39l I4TI 1101] . Consequently, the corresponding charging energy, Ec — e 2 /2C, 
establishes a new relevant energy scale. If the charging energy is larger than the 
thermal energy, Ec 3> k^T ', where T denotes temperature and /cb is the Boltzmann 
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Figure 2. Schematics of ferromagnetic single-electron transistors. The parallel 
and antiparallel magnetic configurations of the system are also specified. The 
island is separated from external electrodes by tunnel barriers. Each barrier is 
characterized by its spin-dependent resistance R ra and capacitance C r (r = L, R). 
The system is symmetrically biased, Vl = V/2, Vr = — V/2, and there is also a 
gate voltage V g applied to the island. 



constant, the effects due to discreteness of charge become observable in transport 
characteristics [T4l 115] , 

The single-electron transistors considered in this section are illustrated 
schematically in Fig. [2j where part (a) shows a device whose all three electrodes 
are ferromagnetic, whereas part (b) shows a system with ferromagnetic source and 
drain electrodes and a nonmagnetic island. Generally, magnetic moments of the 
leads may form an arbitrary magnetic configuration. However, we focuss on systems 
whose moments form either parallel (P) or antiparallel (AP) magnetic configuration, 
as shown in Fig. [2] In the system shown in part (a) the magnetic moments of 
external electrodes point in the same direction, while magnetic moment of the island 
is either parallel or antiparallel to them. In the system shown in part (b), magnetic 
moments of external electrodes in the parallel configuration are aligned, while in the 
antiparallel configuration they are anti-aligned. The two magnetic alignments can 
be easily achieved by sweeping magnetic field through the hysteresis loop, provided 
the respective ferromagnetic components have different coercive fields. One may also 
make use of exchange anisotropy to fix magnetic moment of a particular layer and 
rotate magnetic moment of the second layer with a weak magnetic field. Generally, 
there is a chance that electron tunnelling through a barrier will change its spin 
orientation. However, we consider only spin-conserving tunnelling processes through 
the two barriers. First, we assume that the islands are relatively large, so the effects 
due to quantization of the corresponding energy levels can be neglected. For smaller 
islands, however, the discrete energy spectrum may modify transport characteristics 
[40l 11021 1103| . and this will be considered later (see section [2~3| . 

To describe charge and spin transport we need to specify a model Hamiltonian of 
the system. First, the electrostatic energy required to add n excess electrons to the 
island, while keeping constant voltages Vt and Vr in the left and right electrodes and 
the gate voltage V g , is given by |104j 

H ch =E c (n-Qt) 2 , (5) 
e 

where a constant term (independent of n) is irrelevant and has been dropped. 
Here, the total island's capacitance C is the sum of the capacitances of the left 
and right junctions and of the gate, C = Cl + Cr + C g . The external charge 
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Q g = C'lVl + CrVr + C g Vg accounts for the effect of applied voltages, and can 
be continuously tuned. The total Hamiltonian of the device takes then the general 
form 

H= E H r + H! + H ch +H T = H + H T . (6) 

r=L,R 

Here Hq describes the decoupled leads and island, while Ht takes into account the 
lead-island coupling. The ferromagnetic leads and the island are described by 

Hr ^ ^ ^rkcT^ykCT^Vko" , (7) 
ko- 

for r = L,R, I, where c r k CT are the Fermi operators for electrons with a wavevector k 
and spin a in the electrodes and island (r ~ L,R,I), and e r ko- is the corresponding 
single-particle energy. The last part of the Hamiltonian, 

H T = E E ^kq CT 4 ko .ci qCT + h.c. (8) 

r—L, R kqcr 

describes tunnelling processes between the leads and island, with i,kqo- being the 
relevant matrix elements. 



2.1. Transport in the sequential tunnelling regime 

When the resistances of both tunnel barriers are much larger than the quantum 
resistance, R r ^> Rq = h/e 2 , and the system is not in the Coulomb blockade regime, 
transport is dominated by sequential tunnelling processes. As a consequence, the 
charge is well localized in the island and the orthodox tunnelling theory is applicable 

[mri!]. 

In order to calculate electric current in the sequential transport regime, one may 
use the method based on the master equation which is a detailed balance of electrons 
tunnelling to and off the island. In the stationary case the master equation reads 

o = - E K» + r L» + r £» + r*»] P(n, v) 

a 

+ E K> - !) + r *> - !)] P ( n - !< V ) 

+ E t r L> + !) + r «> + !)] P (" + V ^ ( 9 ) 

where P(n, V) is the probability to find the island in a state with n additional electrons 
when a bias voltage V (V = Vl — Vr) is applied, and r^ tT (n) is the spin-dependent rate 
for tunnelling of electrons with spin a from the lead r to island (upper sign) and from 
the island to lead r (lower sign), when the island is occupied by n excess electrons. 
These tunnelling rates depend on the bias voltage (not indicated explicitly), and can 
be expressed by means of the Fermi golden rule as 

r± r„) = _J AJ ^» ( io) 

rA ' e 2 R ra exp [AE±(n)/k B T] - 1 ' { ' 

where AE^ a (n) describes a change in the electrostatic energy of the system caused by 
the corresponding tunnelling event, when in the initial state there were n additional 
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electrons on the island. In the above equation R ra denotes the spin-dependent tunnel 
resistance of the r-th barrier, given by 

h 

Rra = ~Z 2 U 12' ( ) 

2ne^prcrpia\tr<7\ 

where i r kqo- = t ra has been assumed for simplicity. The spin dependence of the 
resistance is a consequence of the spin-dependent density of electron states at the Fermi 
level in the respective electrodes and the corresponding tunnelling matrix elements. 
In particular, in Eq. (|lip it is due to the spin-dependent density of states of the r-th 
lead, p ra , spin-dependent density of states of the island, pi a , as well as due to the 
spin-dependent tunnelling matrix elements t ra . 

In a stationary state, the net transition rate between the charge states with n and 
n + 1 excess electrons on the island vanishes. The probability P(n, V) can be then 
determined from the following recursion relation J53J HI] 

P(n + 1,V)J2 y a (n + 1) = P(n, V) £ x CT (n) , (12) 

a a 

where x a (n) = Y, r =L,R V tA n ) and y a ( n ) = I2r=L,R T rA n )i corresponding to 
transition rates for tunnelling to and off the island, respectively. 

Generally, energy of an electron after tunnelling event is relaxed to the relevant 
Fermi level in a short time scale. One can assume that the energy relaxation time 
is the shortest time scale, shorter than the time between two successive tunnelling 
events. However, such a restriction cannot be imposed on the spin relaxation time 
which can be relatively long. In a general case, a nonequilibrium magnetic moment 
may accumulate on the island due to the spin dependence of tunnelling processes, 
which leads to spin splitting of the corresponding Fermi level. For arbitrary spin 
relaxation times, the splitting of the Fermi level can be determined from the spin 
balance [M ITO5) 

~{i° R -il)- ei ^^E° F = Q, (13) 

e r sf 

where fij is the island's volume, — e is the electron charge (e > 0), T s f denotes the 
spin relaxation time on the island, while /£ and 1^ are the currents flowing through 
the left and right junctions in the spin channel a. The last term in Eq. (| 13|) takes into 
account intrinsic spin-flip processes on the island. ^From this condition it is possible 
to calculate self-consistently the shifts of the Fermi levels due to spin accumulation 
for both spin orientations. 

The electric current flowing through the left junction can be calculated from the 
following formula: 

oo 

I L =J2 I Z = -eJ2 E [rL(n)-r- La (n)]P(n,V). (14) 

cr (j n— — oo 

Similar formula also holds for Ir. In the stationary state the currents flowing through 
both junctions are equal, It, = Ir = I. 

In the following we discuss two limiting cases: the limit of fast spin relaxation 
and the limit of slow spin relaxation on the island. In the former case the spin of an 
electron tunnelling to the island relaxes before a next tunnelling event takes place. 
In the latter case, on the other hand, the electron spin is conserved for a time much 
longer than the time between successive tunnelling events. 
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Figure 3. Current in the parallel and antiparallel magnetic configurations (a) 
and the resulting TMR (b) as a function of the bias voltage for a FM SET with 
a ferromagnetic island in the limit of fast spin relaxation. The parameters are: 
T = 9 K, = 0.65 MS], = 0.065 MSI, R^ = 5 MSi, 

i?P )V2 (for r = L, R), C L = 0.1 aF, C R = 1 aF, C g 
0, and the offset charge Qo = 0.025e. The dots in part (b) present the 
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experimental data taken from Ref. |39| . 



2.1.1. Fast spin relaxation In the limit of fast spin relaxation, SETs with 
nonmagnetic islands behave like nonmagnetic junctions, exhibiting no TMR effect. 
Therefore, the following discussion in this subsection is limited to FM SETs with 
ferromagnetic islands, see Fig. [2^- Generally, one may assume that the spin-dependent 
resistances for the parallel (-R^o-) an d antiparallel (R^a) configurations fulfill the 
condition R^ a R^ s = R^R^, where a = — a. This formula follows from Eq. (Til) 
when assuming that all spin effects are included into the spin dependent density of 
states (transfer matrix element are independent of spin and magnetic configuration). 
Moreover, since there is no spin accumulation on the island in the limit of fast spin 
relaxation, there is no associated spin splitting of the Fermi level, AEp = 0. 

The current flowing through the system in the parallel and antiparallel 
configurations is displayed in Fig. [3^. For both magnetic configurations, the I — V 
curves reveal the well-known Coulomb steps. Moreover, these two curves are different; 
the current flowing in the parallel configuration is generally larger than the current 
flowing in the antiparallel configuration, see Fig. [3^- This difference leads in turn 
to nonzero TMR effect, as shown in Fig. [3Jd. The TMR effect has a component 
that oscillates as a function of the bias voltage. The amplitude of these oscillations, 
however, decreases as the transport voltage increases. For the parameters assumed 
here, TMR reaches local maxima at the voltages corresponding to the positions of 
Coulomb steps. However, this is not a general rule, and for other parameters TMR 
can have local minima at the Coulomb steps. The global maximum value of TMR in 
Fig. [5b appears at the first step, i.e., at the threshold voltage. When the temperature 
increases, the effects due to discrete charging, i.e. the Coulomb steps and enhancement 
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of TMR at the Coulomb steps become diminished and disappear at k^T ss Eq [98J. 
The oscillatory behavior of the TMR effect with increasing transport voltage was 
observed experimentally, for example by Ernult et al. |39j , see the dots in Fig. [3b- 
The curves presented in Fig. [3] are calculated for the parameters corresponding to 
those in Ref. [39] , To get a good agreement with experimental observations, a nonzero 
offset charge Qq on the island (due to external charges) has been assumed. As can be 
seen in Fig. [3)3, there is a satisfactory agreement between the theoretical curve and 
experimental data for voltages above the threshold, V > 30mV, while in the Coulomb 
blockade regime transport calculated using the sequential tunnelling approximation is 
not properly described. 

2.1.2. Slow spin relaxation In the case of a FM SET with a nonmagnetic island, see 
Fig. [2Jd, a nonzero TMR in the sequential tunnelling regime can exist only when the 
spin relaxation time is sufficiently long, i.e. significantly longer than the time between 
successive tunnelling events. We note that the longest spin relaxation times were 
measured for aluminium and copper [1061 1 107] (for example, the relaxation time for 
copper was estimated to be of the order of 10 -7 s). If this is the case, a nonequilibrium 
magnetic moment builds up on the island due to spin accumulation. This moment 
leads a nonvanishing TMR. In other words, the island becomes magnetized in a 
nonequilibrium situation, and the created moment depends on the bias and gate 
voltages. 

Since the density of states at the Fermi level in a nonmagnetic island is 
independent of the spin orientation, one finds AEp = —AEp. Apart from this, 
the resistances in the antiparallel configuration are = R^ a and — see 
also Fig. Et. 

Typical transport characteristics of a FM SET with nonmagnetic island in the 
limit of slow spin relaxation are shown in Fig. 2] as a function of the bias voltage. The 
splitting of the Fermi level is displayed in Fig. [4^, for both magnetic configurations. As 
one could expect, the splitting takes place only in the antiparallel alignment, while in 
the parallel configuration there is no spin accumulation. In the former case the ratio 
of tunnelling rates for electrons with opposite spin orientations becomes the same 
for electrons tunnelling to and off the island only when a nonequilibrium magnetic 
moment is built on the island. In the parallel configuration, however, this condition 
is already fulfilled without any spin accumulation. However, this is true only in the 
case when both junctions are characterized by equal spin asymmetries. When the 
spin asymmetries of both electrodes are different, spin accumulation also occurs in the 
parallel configuration. 

As illustrated in Fig. [4^, behavior of the Fermi level splitting with increasing 
transport voltage can be decomposed into two components. One component 
monotonously increases, while the second one oscillates with increasing bias voltage. 
This oscillatory behavior can be accounted for in the following way. Let us assume 
that the voltage is slightly above that corresponding to a certain Coulomb step and 
begins to increase. Then, the spin accumulation also increases, until a local maximum 
value is reached. The local maximum occurs at a voltage, at which the chemical 
potential of the depleted spin channel approaches the value which allows the next 
charge state in the island. This, in turn, enhances tunnelling rate (onto the island) 
of electrons corresponding to the depleted spin channel, and consequently reduces the 
spin splitting of the Fermi level. When the voltage increases further, a local minimum 
in the spin accumulation is then reached at a voltage, where the chemical potential 
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Figure 4. Basic characteristics of a FM SET with a nonmagnetic island as a 
function of the bias voltage in the limit of slow spin relaxation on the island: 
(a) spin splitting of the Fermi level, (b) I — V characteristics, and (c) tunnel 
magnetoresistance. The parameters are: k^T = 0.05 Eq, = 5 Mf!, 

Rl i = 2.5 MJi, i?P T = 0.3 Mfi, flg, = 0.15 MQ, C L = C R = C g = 1 aF, 
and V g = 0. 



of the second (accumulated) spin channel approaches the value which allows the next 
charge state on the island. The same scenario repeats at each Coulomb step leading 
to the oscillatory component in the spin accumulation. In the sequential tunnelling 
regime, where only the first-order tunnelling processes are taken into account, spin 
accumulation is exponentially small in the Coulomb blockade regime, as can be seen 
in Fig. |3Jl 

The current as a function of the bias voltage for the parallel and antiparallel 
configurations is shown in Fig.[4}). As before, characteristic Coulomb steps are clearly 
visible. Moreover, owing to different stationary spin accumulations in the parallel and 
antiparallel configurations, the current in the parallel configuration is also different 
from that in the antiparallel configuration. This gives rise to the TMR effect which 
is presented in Fig. Now, the bias dependence is more complex than it was in the 
case of FM SETs with magnetic islands in the absence of spin accumulation. It is 
interesting to note, that TMR can change sign in the vicinity of the Coulomb steps in 
current-voltage curves, which is a consequence of different spin accumulations in the 
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two magnetic configurations. 

Generally, it is more difficult to obey the slow spin relaxation limit in 
ferromagnetic islands than in nonmagnetic ones. Anyway, some experimental data 
show that this is achievable. In such a case, the spin accumulation builds up in 
the magnetic island and has significant influence on transport characteristics. In 
Refs. [93j 105] it has been shown that spin accumulation modifies the 'staircase'-like 
variation of the electric current with the bias voltage. Since the spin accumulation 
depends on the magnetic configuration of the junction, this can also lead to an 
enhanced TMR effect. Moreover, in some voltage regions TMR can change sign. 
In addition, the Coulomb steps for the two magnetic configurations become slightly 
shifted, as was also observed experimentally [49]. The difference between current- 
voltage curves in the fast and slow spin relaxation limits varies continuously with spin 
relaxation rates. This difference was used very recently to evaluate spin relaxation 
time in Co nanoparticles [417] , 

2.2. Contribution to the conductance due to cotunneling processes 

The electric current and spin accumulation in the limit of sequential tunnelling (see 
the preceding subsection) are exponentially small in the Coulomb blockade regime. 
However, charge transport may occur not only due to spin-dependent sequential 
tunnelling, but also due to spin-dependent cotunneling processes [7H [371 1108] , which 
give dominant contribution in the Coulomb blockade regime and also lead to spin 
accumulation [!l08l 1109] , Close to resonance (in the vicinity of the threshold voltage) 
both sequential and cotunneling currents may be comparable. The current I is then 
equal to the sum of first, and second, I^ 2 \ order contributions, I = I 1 - 1 * 1 + I^ 2 \ 

To calculate cotunneling current far from the resonance one could use Eq. 
This formula, however, cannot be used when voltage approaches the threshold voltage 
(at resonance). To calculate current and associated spin accumulation in the whole 
voltage range, the real-time diagrammatic formalism [97j I104[ 1110] 1111] has been 
used. The spin accumulation on the island (or equivalently spin splitting of the 
electrochemical potential) is then determined from the spin balance equation, similarly 
as in the sequential tunnelling regime, see Eq. (|13|) . 

The differential conductance in the whole (sequential and cotunneling) transport 
regime is shown in Fig. [5] for SETs with Ni and Fe electrodes and for nonmagnetic 
islands (in the limit of slow spin relaxation) . Upper part corresponds to the parallel 
configuration while the lower one to the antiparallel one. First, we find well resolved 
splitting of the conductance peaks in the antiparallel alignment, while no splitting 
can be seen in the parallel configuration. This splitting in conductance peaks 
is a direct consequence of the spin splitting of the corresponding electrochemical 
potential of the island, and therefore can be used to detect and measure spin 
accumulation. The absence of conductance splitting in the parallel configuration is 
simply a consequence of the absence of spin accumulation in this configuration. Such 
a splitting of the conductance peaks in the cotunneling regime was recently observed 
experimentally [49] . Generally, there are several experimental techniques by which the 
spin accumulation can be detected indirectly [107[ 1112] . The peculiarities of transport 
characteristics of FM SETs offer new possibilities. 
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Figure 5. Differential conductance in FM SETs with nonmagnetic islands as a 
function of the bias voltage V in the (a) parallel and (b) antiparallel configurations, 
calculated in the whole (sequential and cotunneling) transport regime in the 
limit of slow spin relaxation in the island. The curves have been calculated 
for k^T/Ec = 0.02 and symmetric junctions with Rr| = 5/i/e 2 . The spin 
polarization is p = 0.23 and p = 0.40 for Ni and Fe electrodes, respectively. The 
dotted line corresponds to the sequential tunnelling limit in the SET with Fe 
electrodes. (After Ref. [97]) 



2.3. Discrete energy structure of the island 

Discussion up to now was limited to the case when addition of a single electron costed 
only an electrostatic energy. When, however, the size of metallic island becomes 
reduced further, electron spectrum of the island cannot be considered as a continuous 
one, and discrete structure of the energy levels plays an important role. Addition of 
an electron costs then not only the electrostatic charging energy, but also the energy 
equal to the level spacing AE. In this subsection we describe briefly the effect of 
discreteness on the transport characteristics. The relevant approach was developed in 
Refs. (1021 1103j . The electron transport in the stationary state is then governed by 
the solution of the generalized master equation [1021 1113] 

= - {T(n h n L ) + n Ti ;(n T ,n|) + O ijT (n T , n x )} P(n h n;, V) 

+ T+(n T - 1, n;)P(n T -l,n u V)+ T+(n h nj - l)P(n T , nj, - 1, V) 

+ T-(n T + l,n i )F(n T + l,n i ,^) + r-(n T ,n i + l)P(n T ,n i + l,F) 

+ r! T . x (n T - l,m + l)P(n T - l.n; + 1, V) 

+ n u ( n] - + 1, ni - l)P(n T + l,n x - 1, V), (15) 

where P(n-|-,n;,F) denotes the probability to find n-\ and nj. excess electrons on the 
island (n = n-f + is the total number of excess electrons) . The first term in Eq. (|T5|) 
describes how the probability of a given charge and spin state decays due to electron 
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Figure 6. (Color online) The influence of discrete electronic structure on 
transport characteristics in a FM SET with nonmagnetic island in the limit of slow 
spin relaxation. Voltage dependence of the tunnel current / (a), the differential 
conductance G = dl/dV (b), spin accumulation (n-f — 7i|) (c), standard deviation 
[((rtj — n^) 2 ) — (rif — n l) 2 ] 1 ^ 2 (d)> and tunnel magnetoresistance TMR (e), 
calculated at T = 2.3K. The solid and dashed curves in (a), (c) and (d) correspond 
to the antiparallel and parallel configurations, respectively. The other parameters 
are: AE = 3meV, C L /C L = 5, E c = lOmeV, R L} = 2R Ll = 200MH, .Rrj = 
2R Rl = 4MQ and R m = 2MH for the parallel alignment {2R R ^ = R Rl = 4Mf2 
for the antiparallel alignment). (After Ref. I102p 
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tunnelling to and from the island, whereas other terms describe the rate at which 
this probability increases. The fi-terms account for spin-flip relaxation processes. 
The coefficients entering Eq. (fT5|) are defined as T^(n^,n{) = Yl,r=L R^ra( n T> n l) 
and r(n T ,Ji|) = £ CT [r+ (n T , n L ) + r"(n T , n;)], where r± CT (nt,ru) are the tunnelling 
rates for electrons with spin c, tunnelling to (+) the island from the lead r — L, R or 
backward (— ). These coefficients are given by |102[ 1103] 

r± (nt.ru) =J2^ F -( E i-\ n ^ n l)f ± ( E ^ + E r( n ) - e f), 

% 

n^(n h ni) = ^2^2uj trJ::j9 F+(E la \n h n l )F i r(E jr f\n h n L ). (16) 

i 3 

Here, f + (E) is the Fermi function (/~ = 1 - /+), whereas F+(£i CT |nt, nj (F~ = 
1 — F+) describes the probability that the energy level Ei a is occupied by an 
electron with spin a for a particular configuration (nt,?u). The parameter 
is the bare tunnelling rate of electrons between the lead r and the energy level 
Ei a of the island, and w J(T j"ct is the transition probability from the state ia to ja 
of the island due to the spin-flip processes. The energies E^(n) and E^(n) are 
given by E% (n) = C R /C eV + C/^n) and E%{n) = ~C L /C eV + ^(n) where 
U ± (n) = E c [2(n - n x ) ± 1] and n x = C g V g /e. 

^,From the solution P(n^n^V) of the master equation (Eq. [T5|) . one can 
determine current flowing through the island, 

I L = -eJ2 E P(n h n u V) {ri a (n hni ) -T- La (n hni )} . (17) 

For further discussion we assume that the discrete energy levels Ei a are spin degenerate 
(nonmagnetic situation) and equally separated with the level spacing AE. 

To emphasize the role of spin accumulation let us assume that the intrinsic spin 
relaxation time on the island is long enough to neglect all intrinsic spin-flip processes. 
The corresponding I — V characteristics for the parallel and antiparallel alignments 
are shown in Fig. [5^. In both cases the electric current is blocked below a threshold 
voltage, and a typical 'Coulomb staircase' appears above it, with additional small 
steps due to the discrete levels of the island. The effects due to discrete charging 
and discrete electronic structure are more clearly seen in the differential conductance 
shown in Fig. \6jp, where the small peaks correspond to new discrete levels taking 
part in transport. The difference between the I-V characteristics for the parallel and 
antiparallel configurations is due to a different spin accumulation in both geometries. 
In Fig. [6): we see the average value of the difference between the numbers of spin- 
up and spin-down excess electrons on the island, (M) = (n-f — ni), i.e., the spin 
accumulation. There is no significant spin accumulation in the parallel configuration. 
The number M = n-t — of spins accumulated on the island fluctuates in time 
around its average value (Af), as shown in Fig. [5Ji, where the standard deviation 
((M 2 ) - (Af) 2 ) 1 / 2 = [((n t - ru) 2 ) - (n t - n^ 2 ] 1 / 2 is plotted against the voltage V. 
It is worth to note that although there is almost no spin accumulation in the parallel 
configuration, the corresponding fluctuations are relatively large. 

The difference between the I — V curves in the parallel and antiparallel 
configurations leads to the tunnel magnetoresistance, shown in Fig. [5^. The broad 
peaks correspond to the Coulomb steps, while the fine structure originates from the 
discrete structure of the density of states of the island. 
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One can extend the presented formalism also to the case of a FM SET with 
a ferromagnetic island [114) and calculate nonequilibrium spin accumulation on the 
island (when spin relaxation time is sufficiently long). Using similar approach Inoue 
and Brataas |115j analyzed current-induced magnetization reversal induced by spin 
accumulation rather than by spin torque. They found that the magnetization reversal 
is possible when a free energy change due to nonequilibrium spin accumulation is 
comparable to the anisotropy energy. 

2-4- Shot noise in ferromagnetic SETs 

^From the application point of view, an important transport characteristics of the 
system is the corresponding current noise. Functionality of a system depends on the 
relevant noise to signal ratio, which should be as small as possible. However, the 
noise is also a source of additional information on the system properties like quantum 
and Coulomb correlations [116] . The shot noise in ferromagnetic single-electron 
transistors was studied theoretically in Ref. [117j . where the method developed for 
spinless electrons in single-electron transistors [1181 11191 11201 1121L 1122) was extended 
to magnetic systems. 

The time correlation function of any two quantities X and Y can be expressed 
as [155] 

(X(t)Y(0))= X nlvni P(n' v n[ ] t\n h n l -,0)Y nrtni P°(n h n i ), (18) 

rij ;rif ,n| 

where, P(n'p n',; t\n-\, nj,; 0) is the conditional probability to find the system in the 
final state with n'^ and n\ excess electrons at time t, if there was n-\ and nj, excess 
electrons in the initial time t=0. The probability P° can be determined from Eq. (j!5[) . 

Following this procedure, the current-current correlation function in a FM SET 
considered in the preceding subsection has been calculated in Ref. [117j . The 
corresponding Fourier transform can be presented as 

S n (u) = S? I h + Sfr{u), (19) 

where SfJ 1 is the Schottky value (the frequency independent part), while the second 
term in Eq. ([19p is the frequency dependent component. 

Figure [7^ shows the bias dependence of the zero- frequency current noise Sjj{cj — 
0). The current noise is smaller in the antiparallel configuration than in the parallel 
one. This is because in the presence of spin accumulation (which is significant only 
in the antiparallel alignment) the amplitude of fluctuations is smaller. In Fig. [7)3 
Sn(u) = 0) is split into two parts; the frequency independent component Sfj 1 and 
the contribution Sfj (u> = 0) arising from the frequency dependent part of the current 
noise, see Eq. (fTT))) . The component Sfj 1 is almost constant w 2eI(C\ 4- C|)/C 2 at the 
plateaux of the I — V curve and increases with opening of new channels. Dynamical 
correlations between the currents are described by Sjj(uj). Its value in the limit uj — > 
can be positive between the I—V steps and negative when new channels become open. 
This is evident for the antiparallel alignment at V ~ 26mV, when opening a tunnelling 
channel for electrons with a =J, leads to negative dynamical correlations. This effect 
is almost compensated by an increase in Sj/\ and therefore one gets only a small 
reduction of the current noise Sn (u> = 0). 

In the power spectrum of the current Sjj(u>) one can distinguish two distinct 
relaxation times, one in the high and another one in the low frequency regions [117] . In 
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Figure 7. Voltage dependence of the current shot noise at ui = (a) in a 
FM SET with nonmagnetic island and in the slow spin relaxation limit for 
4H LT = R L i = 8MC, R R ^ = 4R Ri = 240MO for the antiparallel alignment 
(iii-t = &Rli = 8MS7, for the parallel configuration) (other parameters as in 
Fig-|6]l. In part (b) Sii(u = 0) is split into two components: Sfp (upper curves) 
and SfJv = 0) (lower curves). (After Ref. [TTZ]) 



a wide voltage range the corresponding relaxation times are very close to the effective 
relaxation times for the charge and spin noise. 

The asymmetry between the tunnelling channels for electrons with the opposite 
spins leads to activation of the spin component in the current noise. In Ref. [117 
the components Sjj charge and Sfj spin , corresponding to the charge and spin noise, 
respectively, have been extracted from Sfj. It has been shown that the charge 
component is almost constant whereas the spin component increases with spin 
polarization p of the leads, and for p — * 1 can be much larger than the charge 
component. The analysis showed that both charge and spin fluctuations are relevant 
for the shot noise in FM SETs. It has been also pointed out that a super-Poissonian 
shot noise, the Fano factor Sjj(Q)/2eI > 1, can occur due to lifting of the spin 
degeneracy. 

3. Transport in double-island devices 

Some new features of transport characteristics appear in spin-polarized electronic 
transport through double-island structures [351 S01 [Ml (Ml I124J . A typical double- 
island device, shown schematically in Fig. [51 consists of two metallic islands separated 
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Figure 8. Schematic of a double-island device. The first (1) and second (2) 
islands are capacitively coupled to two gate voltages, V g \ and V g 2, and separated 
from each other and from the left and right electrodes by tunnel barriers. 



from each other and from external reservoirs by tunnel barriers. A gate voltage 
is additionally attached capacitively to each island. Generally, both islands and 
external electrodes can be either magnetic or nonmagnetic [UJ I3H1 H25| . In the 
following, we consider the situations when at least two electrodes (external or central) 
are ferromagnetic, and their magnetic moments form only collinear (parallel and 
antiparallel) configurations. Apart from this, we limit discussion to continuous density 
of states in the islands. 

Hamiltonian of the double-island systems is similar to that used for single-island 
devices, but now the electrostatic energy is given by [1241 1126] 

9A\^ f n ,-9^Y 



H c h = E Cl yni —J + E C2 yn 2 ^ ( 

+ 2E Cm (m - [n 2 - , (20) 

where is the number of excess electrons on the first (second) island and Q 9 i( S 2) 
describes the respective charge induced by applied voltages, Q g i( g 2) = Cl(r)Vl(r) + 
C g i(g2)V gl ( g 2)- Apart from this, Ec x and Ec 2 denote the charging energies of the 
corresponding islands, whereas Ec M describes electrostatic coupling of the islands, 



E C 1 (C 2 ) - 7^ — 1- 7T7r ' ( 21 ) 
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Ec M = wr - h^- 1 . ( 22 ) 

with C\(2) being the total capacitance of the first (second) island, Cu2) — Ci(ii) + 
Cgi( g 2) + Cm, and Cm denoting the capacitance of the middle junction. 

To find the stationary current flowing through the system one needs then to know 
the probabilities P(rii, ri2, V) of finding the system in the charge state with n\ and n 2 
additional electrons on the first and second islands, respectively, when a bias voltage V 
is applied. These probabilities can be calculated in a recursive way from the following 
steady-state master equation 

o=-E[riiK« 2 ) + r^(m,n 2 ) + r^(m,n 2 ) 

a 

+ r 2 r 1 (ni,n 2 ) +r 2 r R (ni,n 2 ) + T^C^i) "2)] n 2 , V) 
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+ r "( ni _ - 1, n 2 , V) 

a 

+ J2 r ii( ni + X ' + 1, 712, V) 

(7 

+ J2 r i2( n i + i. »*2 - + i. «2 - 1, V) 

+ 51 r 2i(™i - !> «a + - 1, n a + 1, V) 

(7 

+ E r- 2 (m, n 2 - l)P(ni, n 2 - 1, V) 

a 

+ ^2 T 2 R (ni,n 2 + l)P(n 1 ,n 2 + l,V), (23) 

a 

with the normalization condition ^ n „ 2 P{rii, n 2 , V) = 1. The corresponding 
transition rates are given by Eq. (|10|) , with the respective changes in the total system 
electrostatic energy. For example for tunnelling from the first island to the second 
one, the change in the electrostatic energy can be written as 

AE^ 2 (m,n 2 ) = E(m - l,ri2 + 1) - E{n x ,n 2 ) - AE F1 + AE F2 . (24) 

Here, AEp 1 and AE F2 denote the corresponding shifts of the chemical potentials for 
spin a in the first and second islands, respectively. The spin asymmetry of tunnelling 
through the barrier r will be characterized in this section by the asymmetry factor 
a r , defined as a r = R r ^/R r ^. This asymmetry factor corresponds to the ratio of 
the respective spin-dependent densities of states. In particular, for FM/NM junctions 
it is given by a = Pi/p^, whereas for FM/FM junctions a — p 2 ^/pp provided that 
the two ferromagnetic electrodes are built of the same material. We note that the 
relation between spin dependent barrier resistances in the parallel and antiparallel 
configurations of magnetic moments on the opposite sides of the barrier is the same as 
that described in Sec. 2. When a bias voltage is applied to the system, a nonequilibrium 
spin accumulation may appear on the islands. Generally, the shifts of the Fermi 
level for both spin orientations are different. However, one can assume that the 
ratio of the Fermi level shifts for the spin-up and spin-down electrons, defined as 
/3j = —AE^pj/AEpj, fulfills the relation (3j = piji/piji, for the first (j — 1) and 
second (j = 2) island, respectively, with pij a being the spin-dependent density of 
states of the j-th island. As a consequence, for nonmagnetic islands one directly gets, 
p 1= P 2 = l. 

Electric current flowing through the system can be determined from the following 
formula: 

oo 

lL=-eJ2 J2 [rL(m,n 2 )-r^(ni,n 2 )]P(n 1 ,n 2 ,y), (25) 

a ni,U2 — — oo 

which corresponds to the current flowing through the left junction, but 1^=1^ = 1 
in the stationary state. 

3.1. Fast spin relaxation: no spin accumulation 

Consider first transport characteristics of a system built of two ferromagnetic 
islands and nonmagnetic external electrodes, shown in Fig. [3J Different magnetic 
configurations of the system are specified in the inset of Fig. [9};. When the spin 
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Figure 9. The average electron number on the islands (a), currents (b) in 
the parallel and antiparallel configurations and the resulting TMR (c) as a 
function of the bias voltage. The parameters are: Ci = Ci = C = 3Cm = 3 
aF, k^T/(e 2 /2C) = 0.05. The spin asymmetries of resistances in the parallel 
configuration are aj, = an = 5, «m = 25, whereas the total junction resistances 
are ij£ = R R = R M /10 = 1 Mf2. In the antiparallel configuration, R^ = fi£ CT , 
pAP _ nAP _ CrP dP U/2 ,_j pAP _ nP 



relaxation time in the islands is much shorter than the time between two successive 
tunnelling events, no spin accumulation builds up on the islands. Since there is an 
asymmetry between the two barrier resistances, some charge accumulates on the 
islands, as displayed in Fig. [9^- First of all, the magnitude of excess charge on 
the islands increases with increasing voltage. For the barrier asymmetry assumed 
in Fig. [9j the electrons easier tunnel to the second island from the right lead than 
out of the second island to the first one. As a consequence, the electrons accumulate 
on the second island. On the other hand, the electrons easier tunnel out of the 
first island to the left lead than from the second island to the first one. Thus, 
the number of excess electrons on the first island decreases with increasing the bias 
voltage - there are holes accumulated on the first island. This occurs in both magnetic 
configurations. The electric current flowing through the system in the parallel and 
antiparallel configurations is shown in Fig. [SJa. The Coulomb steps due to discrete 
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Figure 10. The shifts of the Fermi levels for spin-up electrons (a), currents in the 
parallel and antiparallel configurations (b) and TMR (c) as a function of the bias 
voltage. The parameters are: C\ = C 2 = C = 3C M = 3 aF, k B T/(e 2 /2C) = 0.05, 
01 = (h = 1, r s /i — » oo, T s f 2 =—* co, a L = a R = 5, a M = 1, whereas 
RP = i?P = R^/50 = 1 Mf2. In the antiparallel configuration, = R^ a , 

dAP _ pP 



charging are clearly evident, and the difference in currents flowing in the parallel and 
antiparallel configurations gives rise to the tunnel magnetoresistance, shown in Fig. [5J;. 
The TMR effect oscillates as a function of the bias voltage and the amplitude of these 
oscillations decreases as the voltage is increased in a similar way as in the case of 
single-island FM SETs, see Fig. [3)3. Moreover, some dips occur now in TMR at the 
voltages corresponding to the steps in the current-voltage characteristics. 

3.2. Slow spin relaxation: spin accumulation 

When the spin relaxation time is longer than the time between two successive 
tunnelling events, a nonequilibrium magnetic moment appears on each island. The 
corresponding shifts of the Fermi level due to spin accumulation can be calculated 
from spin current conservation, see Eq. (| 13[) . written for each island |124j . 

Figure [10] shows the shifts of the Fermi levels and currents in the parallel and 
antiparallel configurations, as well as the resulting TMR calculated in the limit of long 
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Figure 11. The shifts of the Fermi levels for majority-spin electrons (a), currents 
(b) and differential conductance (c) in the parallel and antiparallel configurations, 
and the resulting TMR (d) as a function of the bias voltage. The parameters are: 
T = 140 K, C L = 0.45 aF, C M = 0.2 aF, C R = 0.35 aF, C gl = C g2 = 0, 
a 2 L = olm = olr = 25, and Pi = P2 = 0.2, whereas T 3 fi = r 3 f2 = 00. The total 
junction resistances are Rl = 3500 MQ, Rm = Rr = 1 MO. In the antiparallel 
configuration, R^ p = R? , and iJ A t p = R A f = (R^RF.) 1 / 2 , for r = M,R. 
(After Ref. [T24] 1 



spin relaxation time for the system built of ferromagnetic electrodes and nonmagnetic 
islands. The magnetic moments of external electrodes can form either parallel or 
antiparallel configurations, as illustrated in the inset of Fig. [TUb . First of all, the 
nonequilibrium spin accumulation, shown in Fig. llOh . exists not only in the antiparallel 
configuration but also in the parallel one. In the antiparallel configuration the shifts 
of the Fermi level for a given spin orientation are equal on both islands, whereas in 
the parallel configuration they are opposite. The effects due to discrete charging lead 
to an oscillatory behavior of the Fermi level shift, in a similar way as in the case 
of single-island FM SETs with nonmagnetic islands discussed above. The currents 
flowing through the system in both magnetic configurations are shown in Fig. 110b . 
Due to nonequilibrium spin accumulation induced on the islands, these currents are 
different, which leads to nonzero tunnel magnetoresistance, as displayed in Fig.flOb. It 
is interesting to note, that now TMR changes sign in certain transport voltage regions. 
These effects are clearly due to magnetic moments accumulated on the islands. If the 
spin relaxation time becomes shorter than the time between successive tunnelling 
events, spin accumulation disappears and, consequently, TMR also vanishes. 

Recently, several experiments on spin-polarized transport through granular 
systems were reported [33 [Ml US HH W\- For example, in Ref. [52 , the 
tunnelling current was driven from a tip of scanning tunnelling microscope through 
ferromagnetic grains to ferromagnetic electrode. Transport measurements of such 
devices showed pronounced Coulomb steps in the current-voltage characteristics. 
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Figure 12. (Color online) Differential conductance in the antiparallel 
configuration as a function of the bias voltage calculated for different values of the 
spin relaxation time T s /i = T s f2 = r B f and for pf-^Qji = pf 2 ^l2 = 1000/eV. The 
other parameters as well as magnetic configuration of the system are the same as 
in Fig. [HI (After Ref. [124]) 



Moreover, negative differential conductance was observed. It was further proposed 
that the negative differential conductance could be a consequence of a nonequilibrium 
spin accumulation. Such systems can be modelled theoretically by double-island 
devices whose two islands and the right electrode are ferromagnetic, whereas the left 
electrode is nonmagnetic, corresponding to the nonmagnetic tip of scanning tunnelling 
microscope, see the inset of Fig. [TTh . Transport characteristics of such a device are 
displayed in Fig. [TTJ for the parameters taken from Ref. [92] . The shifts of the Fermi 
level for spin-up electrons are shown in Fig. [TTa . These shifts arc different in both 
magnetic configurations. Generally, spin accumulation on the first island is larger 
than the accumulation on the second island. The reason for this is the fact that 
the rate for electron tunnelling from the first island to the left lead is smaller than 
the rate for tunnelling of electrons to or from the second island, which is due to 
asymmetry of barriers. The currents flowing through the system in the parallel and 
antiparallel configurations are illustrated in Fig. II lb . Moreover, negative differential 
conductance occurs in both magnetic configurations, however, it is more pronounced 
in the antiparallel configuration, as shown in Fig. [1Tb . It can be seen that negative 
differential conductance increases with increasing the bias voltage. The resulting TMR 
effect is displayed in Fig. [TTH . It is interesting to note that TMR oscillates between 
negative and positive values. 

Oscillations of the sign of TMR and differential conductance result from spin 
accumulations in both islands and are absent in the limit of fast spin relaxation (no 
spin accumulation) in the islands. The results presented in Fig. [TTJ correspond to the 
long spin relaxation limit. In such a limit, some spin accumulation may occur even for 
a very small current flowing through the system, giving rise to NDC and oscillations of 
the TMR sign for small bias voltages. However, both effects disappear when the spin 
relaxation time is shorter than the time between successive tunnelling events. Thus, 
for a finite relaxation time, one may expect absence of NDC and TMR oscillations 
for small voltages and the onset of these effects at larger voltages. This is because at 
some voltage there is a crossover from the fast to slow spin relaxation limits. In fact 
such a behavior is consistent with experimental data of Ref. [92] . 
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Figure 13. (Color online) TMR as a function of the bias voltage calculated 
for different values of spin relaxation time T s fi = T B f2 = r B f and for p^Sl/i = 
pf 2 flj2 = 1000/eV. The other parameters as well as magnetic configuration of 
the system are the same as in Fig. 1111 (After Ref. |124| ~l 



The disappearance of NDC with decreasing spin relaxation time r s f is shown 
explicitly in Fig. [12l where the bias dependence of differential conductance is presented 
for different values of the spin relaxation time. This figure clearly shows that 
NDC disappears when spin relaxation time decreases, in agreement with the above 
discussion. In the limit of fast spin relaxation time, the differential conductance is 
positive, although its periodic modulation still remains. 

Similarly, periodic oscillations of the sign of TMR also disappear with decreasing 
spin relaxation time r s f . This behavior is shown in Fig. 1131 where the bias dependence 
of TMR is shown for several values of r s f. First, the transitions to negative TMR 
disappear with decreasing r s f. The TMR becomes then positive, although some 
periodic modulations survive. Second, the phase of the modulations shifts by about 
7r when the spin relaxation varies from fast to slow limits. 

In the case of analyzed systems, the negative differential conductance occurs due 
to nonequilibrium spin accumulation on the islands. It is however worth noting that 
the negative differential conductance may also exist in single-electron devices built of 
nonmagnetic materials |127| 1128] . 

4. Spin polarized transport through single-level quantum dots connected 
to ferromagnetic leads 

In the previous two sections we considered spin-dependent transport through single- 
electron devices, whose central electrodes (islands) were described by continuous 
(or discrete with small level separation) energy spectrum, and the most relevant 
and dominant energy scale was the electrostatic charging energy. In ultra-small 
metallic islands (nanoparticles) or in semiconducting quantum dots, the level spacing 
is comparable with the charging energy or even larger. This also happens in the case 
of a molecule attached to metallic leads. In this limit one arrives at slightly more 
sophisticated single-electron devices jTH I129[ I130j . 

In this section we consider a FM SET based on a semiconductor quantum dot 
coupled to ferromagnetic leads. This model also applies to molecules attached to 
ferromagnetic electrodes. Single-electron transistors based on quantum dots are 
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of current interest not only because of new and interesting physics emerging in 
those systems, but, more importantly, due to possible future applications and due 
to the possibility of manipulation of a single electron charge and a single electron 
spin [T51 33 11311 1132] , Furthermore, quantum dots are also interesting for future 
applications in quantum computing |133[ 1134) . 

Transport properties of quantum dots coupled to nonmagnetic leads have already 
been extensively studied both theoretically and experimentally [HI [66l Ell E21 
I131| I132|, I135|, 1136] . However, further interesting effects occur in the case of 
quantum dots coupled to ferromagnetic leads, e.g. spin accumulation, parity effect 
on tunnel magnetoresistance, zero-bias anomaly in the Coulomb blockade regime, 
exchange field, splitting of the Kondo anomaly, and others. Most of the works 
concerned theoretical description of spin-polarized transport in the weak coupling 
regime, as well as in the strong coupling regime, where the Kondo physics emerges 



[1371 11351 EH EDI EH HH HH HH [US] . Sequential transport through a single-level 



quantum dot coupled to ferromagnetic leads was studied for both collinear [1461 1147] 
and non-collinear [1481 11491 11501 I151[ I152L I153j configurations of the electrodes' 
magnetic moments. Spin-polarized transport in the cotunneling regime has also been 
addressed for collinear systems |154[ I155| I156| 1157] , as well as for systems magnetized 
non-collinearly [158[ 1159"] 1160} fTHTj . Furthermore, the resonant tunnelling was also 
considered [1321 PTES] . 

Quantum dots coupled to ferromagnetic leads may be realized experimentally 
in various ways, including ultrasmall metallic (e.g. aluminum) nanoparticles [41j, 
single molecules [51], granular structures |163j . self-assembled dots in ferromagnetic 
semiconductors [46] , carbon nanotubes [52 ] [53 ] [54 ] [55 ] [56 ] [57 ] 158 ] 159 ] l60l I6T] [62] . and 
magnetic tunnel junctions [39]. Quite recently, semiconductor quantum dots based 
on two-dimensional electron gas were successfully attached to ferromagnetic leads 



4-1- Quantum dots weakly coupled to ferromagnetic leads: collinear magnetizations 

Let us begin our discussion with the case when the dot is weakly coupled to the leads. 
One can then use a perturbative approach in which the coupling is considered as a 
small perturbation to the unperturbed quantum dot and the leads. Transport is then 
dominated either by the first order or by the second order (cotunneling) contributions. 
Later on we consider a more general situation when the coupling may be strong and 
a new Kondo physics emerges at low temperatures. 

4- 1.1. Model and method A schematic of a quantum dot coupled to ferromagnetic 
leads is presented in Fig. [T3] The magnetizations of the leads can be either parallel or 
antiparallel. There is also a gate voltage attached to the dot. The system is modelled 
by an Anderson-like Hamiltonian of the general form [164] 



where the first and second terms describe the left and right reservoirs of noninteracting 
electrons, H r = ^ £rq<TcJ q<T c rq(T for r = L, R [see Eq. ([7])]. The third term of the 
Hamiltonian, H^, represents the dot and includes two components: one describes 
noninteracting electrons in the dot level and the other represents the Coulomb 
interaction of two electrons residing in this level, 



nana mug. 
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Figure 14. Schematic of a quantum dot coupled to ferromagnetic leads. The 
magnetic moments of external electrodes can be aligned either in parallel or 
antiparallel. The system is symmetrically biased and there is a gate voltage 
attached to the dot. 



where U is the correlation energy, e a — e =F A/2 is the energy of an electron in the dot 
with spin a, and dj. (d a ) is the corresponding creation (annihilation) operator. The 
position of the dot level can be tuned by the gate voltage, but is independent of the 
symmetrically applied transport voltage. Generally, the dot level may be spin-split, 
for example due to a stray field of the electrodes or due to an external magnetic field. 
The corresponding level splitting is denoted by A, whereas e is the energy of the spin- 
degenerate dot level. In general, four different states of the dot are possible: empty 
dot (x = 0), singly occupied dot with a spin-up (\ =T) or spin-down (x =4) electron, 
and doubly occupied dot (x = d), where |x) are the corresponding eigenfunctions. 

Interaction between the leads and quantum dot is incorporated in the tunnelling 
Hamiltonian, Ht, given by 

H T = (W4q<A + KqadUr^) , (28) 

r=L,R qer 

where i rqCT are the tunnel matrix elements. Tunnelling gives rise to an intrinsic 
broadening T a of the dot levels, T a = J2 r =LR^r- The parameters rj and r| 
describe contributions to the level widths due to coupling of the dot to the lead r. 
The respective contribution can be expressed in terms of the Fermi golden rule 
as FJT = 27r^ q |£ rqcr | 2 <5(w — £ r qer)- Assuming the tunnel matrix elements t rcvJ to be 
independent of the wave vector q, one can write 

Y° = 2T<\t ra \ 2 Pra . (29) 

The coupling parameters are usually expressed in terms of the spin polarization p r 
of the lead r, defined by Eq. ©, as rl {i) = T r {l ± p r ), where T r = (rj + T l r )/2. 
As reported in [165] . typical values of the dot-lead coupling strength T in the weak 
coupling regime are of the order of tens of /xeV. 

In order to investigate transport properties of the system in the whole range of 
parameters, one has to use a method which is more sophisticated than that based 
on the Fermi golden rule. Two such techniques are commonly used: the method 
based on the equation of motion for the electron Green functions, and the real-time 
diagrammatic technique. Discussion in this section is based mainly on the latter one 
[1041 11661 11671 11681 11691 1170j . This technique is based on a systematic perturbation 
expansion of expectation value of the current operator and the density matrix order 
by order in the dot-lead coupling strength T. Even, if we limit further considerations 
based on this technique to the second order processes, this allows us to describe 
transport in the resonance regime, where the usual second order perturbation term 
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diverges. One also could go beyond the second order theory and describe more subtle 
effects. 

Within the real-time diagrammatic technique, the density matrix elements of the 
quantum dot, P$£, are given by a kinetic equation in the Liouville space 

° = fe«-e*.)^ + £*£^. (30) 



X2X2 X2 

xix 2 



where X } is the irreducible self-energy corresponding to transition forward in time 

from state \x\) to \xi) and then backward in time from state \xs) to | X2) - 

When the tunnelling processes are spin-conserving and magnetic moments of the 
leads are collinear, the density matrix is diagonal. After performing the perturbation 
expansion, one gets the following first-order (sequential tunnelling) and second-order 
(cotunneling) master equations: 

o=£ s M 0) > (3i) 

^EW+W. ( 32 ) 

x 

respectively, where the probabilities obey the normalization condition P^ 71 ^ — 

Systematic perturbation expansion of electric current with respect to the coupling 
strength T can be performed in a similar way, and the first-order and second-order 
contributions to current are given by the expressions [1661 1169] 

'^-IMM^ (33) 

xx' 

j( 2 ) = _!ly- s i(2) P (o) + s i(i) P (i) (34) 

2ii X 'X X x ^x'x X ' 

xx' 

where the coefficients S^'x an< ^ ^x'x are nrs ^" an< ^ second-order self-energies, 
modified as compared to ^v^, to account for the number of electrons transferred 
through the barriers [156] . 

Transport characteristics of quantum dots are generally studied in both the linear 
and nonlinear response regimes [156] . It is thus important to distinguish between 
different transport regimes of quantum dots, which are sketched in Fig.[T51and labelled 
by the corresponding capital letters. 

First of all, by changing position of the dot level (by the gate voltage for instance) 
or applying the bias voltage, one can cross over from one regime to another. The three 
regions (A, B and A') around zero bias correspond to the regime where sequential 
tunnelling is exponentially suppressed and current flows mainly due to cotunneling. 
The charge state of the dot is then fixed (strictly at zero temperature) to zero electrons 
in the regime A, one electron in the regime B, and two electrons in the regime A'. The 
first-order tunnelling processes are possible once the bias voltage is increased above 
the threshold voltage, allowing for finite occupation of two adjacent charge states (zero 
and one for regime C, and one and two for regime C). In the regime D all four dot 
states x = 0,T,l,d are possible. By performing a particle-hole transformation, the 
behavior in regime A' and C can be mapped to that in regime A and C, respectively. 
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Figure 15. A sketch illustrating different transport regimes. The respective 
regimes are separated by solid lines and labelled correspondingly. 



4- 1.2. Case of nonmagnetic leads To discuss the role of second-order processes 
we first briefly describe the case of a quantum dot coupled to nonmagnetic leads 
(pl = PR = 0). Figure [TBI shows the first-order (dashed line) and second-order (dotted 
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Figure 16. Linear conductance as a function of the level position. The dashed 
line corresponds to the first-order contribution G", the dotted line presents the 
second-order conductance and the solid line shows the sum GW +G^ 2 \ The 
different transport regimes are also specified. The parameters are: k^T = T, 
U = 20r, and p = 0. (After Ref. [156]) 

line) contributions to the linear conductance as well as the total conductance (solid 
line). The conductance is shown there as a function of the dot's level energy. When 
the dot level crosses the Fermi level of electrodes, there is a resonance peak in the 
linear conductance. Another resonance appears when e + U crosses the Fermi level. 
The resonance peaks acquire a certain width as a result of the level broadening due to 
coupling to the leads (thermal fluctuations also contribute). It is interesting to note, 
that the second-order contribution becomes negative at resonances, which indicates 
that the second-order processes renormalize the first-order (sequential) contribution. 
Except for resonances, the dot is either in the empty (regime A) or doubly occupied 
(regime A') state, or in the Coulomb blockade regime (regime B). In all these three 
cases the cotunneling contribution to electric current becomes dominant. It is also 
worth noting, that the second-order processes lead to renormalization of the dot level 
energy [166] . 
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4-1.3. Case of ferromagnetic leads When the external electrodes are ferromagnetic, 
qualitatively new features appear in the transport characteristics. These are 
particularly pronounced in the tunnel magnetoresistance, defined by Eq. |T]), Below 
we describe some of them. 
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Figure 17. (Color online) The first-order (a) and first- plus second-order (b) 
tunnel magnetoresistance as a function of the bias and gate voltages for the 
parameters: k B T = 1.517, U = 40I\ and p = 0.5. (After Ref. [T561 1 



Tunnel magnetoresistance calculated in the first order of the dot-lead coupling 
strength T (sequential transport limit) is shown in Fig. [T"7r as a function of the bias 
and gate voltages. It is evident that the TMR generally acquires then two different 
values, depending on the transport region. For the regions A (and A'), B and D, the 
TMR value is 

TMR*' B,D = = iTMR Ju11 , (35) 

H f — p l 2 

while for the region C (and C) it is 

TMR s c eq = 3(i V p2) = ^TMR Ju11 . (36) 

This behavior can be accounted for as follows. For the regions A (A') and B the 
first-order linear conductance in the parallel and antiparallel configurations is 

Gp 13 ~ r/2 and G^ P ~ T(l-p 2 )/2 , (37) 

which leads to the magnetoresistance equal to TMR = TMR^ c ' q B = p 2 /(l — p 2 ). To 
account for the behavior of TMR in the regions C and D let us consider the zero 
temperature limit. In the region C there are then three dot states taking part in 
transport: x = 0>1\l (because of the particle-hole symmetry, the results are also 
applicable to regime C). On finds then the first-order currents in the parallel and 
antiparallel configurations to be J p ~ T/3 and J^p ~ r(l— p 2 ) / (3+p 2 ) , which leads to 
TMR s c oq = 4p 2 /3(l-p 2 ). Similar analysis for the region D gives TMRf cq = p 2 /(l-p 2 ). 

As we already know from previous sections, first-order (sequential) transport does 
not describe properly the blockade regions (the regions A, A' and B). Transport in 
these regions is expected to be strongly modified by the second-order (cotunneling) 
processes. Accordingly, the above picture is also strongly modified. The corresponding 
TMR is displayed in Fig. [T7b . Since the sequential tunnelling dominates transport 
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above the threshold voltage (outside the blockade regions, i.e. in the regions C and D), 
the total TMR is only slightly modified there as compared to that in the first-order. In 
the regions A, A' and B, however, the first-order tunnelling processes are exponentially 
suppressed and the current is dominated by the second-order (cotunneling) processes. 
As follows from Fig. [T7b . cotunneling has then a significant influence on TMR. Apart 
from this, TMR in the regions A and A' behaves differently from that in the region B. 
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Figure 18. The total linear tunnel magnetoresistance (solid line) as a 
function of the level position. The dashed line represents the first-order tunnel 
magnetoresistance, whereas the dotted line corresponds to the Julliere's value. 
The parameters are k B T = 1.5T, U = 40r, and p = 0.5. (After Ref. [156]) 

Let us look now in more details at the vertical cross-section of Fig. [T7b along 
the zero bias (linear response). The corresponding linear TMR as a function of the 
level position (gate voltage) is shown in Fig. [TH By changing position of the dot 
level, one crosses over from the regime A through the regime B to the regime A'. For 
comparison, the Julliere's value of TMR, as well as the first-order contribution to TMR 
(denoted as TMR*- 1 - 1 ) are shown there. Note, that the first-order TMR is independent 
of the gate voltage and is equal to a half of the Julliere's value. First of all, one can 
note that the second-order processes modify the TMR substantially. Unlike the first- 
order term, the total TMR (first-order plus second-order contributions) does depend 
on the level position. The interesting feature of TMR shown in Fig. [T5b is a strong 
parity effect. The TMR reaches maximum when there is an even number of electrons 
on the dot (zero for the regime A or two for the regime A'), and has minimum for 
an odd (one in the regime B) number of electrons. A universal feature is that for 
even electron number, TMR exactly coincides with the Julliere's value. The system 
behaves then like a single ferromagnetic tunnel junction. Such a situation can take 
place when tunnelling processes are coherent. The only second-order processes that 
contribute to conductance in these regions are the non-spin-flip cotunneling processes, 
in which the electron spin is conserved. Such processes indeed are fully coherent. The 
corresponding cotunneling rates are proportional to the product of the density of states 
of the left and right leads, thus, one can express the second-order linear conductance 
as 

4 2) ~y(l+P 2 ) and Gi 2 p ~^(l-/), (38) 
for the parallel and antiparallel configuration, respectively. As a consequence, the 
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TMR is then equal to that of a single planar ferromagnetic tunnel junction 

TMR A = Y^p. = ™ rJuU • ( 39 ) 

Because of the particle- hole symmetry, the same result can be obtained for regime A'. 

The situation becomes, however, more complex for an odd number of electrons 
on the dot. Apart from the non-spin-flip processes, there are also spin-flip cotunneling 
processes that change spin of the dot. Such second-order processes give rise to spin 
relaxation on the dot and lead to reduction of TMR in the regime B. The dependence 
of TMR on the position of the energy level in the regime B reflects the relative relation 
of the spin-flip to non-spin-flip cotunneling processes. The minimum value of TMR 
appears for e = —U/2, see Fig. [T8b. where one finds 

™<in = 3(1 2p2 2) = ^TMR JuU . (40) 
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Figure 19. The total current (a) in the parallel (solid line) and antiparallel 
(dashed line) magnetic configurations as a function of the bias voltage in units of 
Io = eT/H. Part (b) shows the first-order contribution to the TMR (dashed line) 
and the total TMR (solid line). The parameters are: k B T = 1.5T, e = -20F, 
U = 40r, and p = 0.5. (After Ref. [156]) 

Electric current flowing through the system in the nonlinear response regime is 
shown in Fig. 119b for both magnetic configurations and for a symmetric Anderson 
model (e = —U/2). The corresponding TMR is also shown there, see Fig. [T^b . 
The contribution due to the second-order tunnelling processes is significant (or even 
dominant) in the Coulomb blockade regime. The sequential TMR is shown by a 
dashed line in Fig. [T^b . One can see that sequential TMR is constant as a function 
of the bias voltage. However, this is not a universal behavior and occurs only for 
symmetric Anderson model. Generally, first-order TMR depends on the applied bias 
voltage, as illustrated in Fig [T7k . The second-order processes lead to a strong and 
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Figure 20. (Color online) The total tunnel magnetoresistancc as a function of 
the bias and gate voltages in the presence of external magnetic field for (a) A = 4P 
and (b) A = -4I\ The parameters are: k B T = 1.5V, U = 40I\ and p = 0.5. 



nontrivial dependence of the total TMR on the bias voltage. For large values of 
the transport voltage (regime D), the first-order tunnelling processes dominate, and 
therefore the total TMR is only slightly modified as compared to that in the sequential 
transport regime. However, for voltages below the threshold voltage (regime B), 
the cotunneling processes lead to a strong bias dependence of TMR in the Coulomb 
blockade regime, and to a deep minimum of TMR in the zero-bias limit. In the case 
of metallic islands, cotunneling processes usually lead to an enhancement of the TMR 
effect in the Coulomb blockade regime [74]. Here, we have the opposite situation, i.e., 
suppression of the effect. At low bias voltage, \eV\ <C feT, the single-barrier spin-flip 
processes reduce the TMR. This is however no longer the case for nonlinear response 
regime, \eV\ 3> kBT, where the spin accumulation diminishes the amount of spin-flip 
processes and the TMR increases. Consequently, the TMR effect in regime B increases 
with rising the bias voltage within the limits TMR Jull /3 < TMR B < TMR Ju11 . The 
minimal value is reached at V = and e = —U/2, whereas the maximal value is 
approached for bias voltages large as compared to thermal energy but still far away 
from the onset of sequential tunnelling. 

New features of transport characteristics appear when the spin degeneration of 
the dot level is lifted, s\ ^ ej,, e.g. due to an external magnetic field applied to 
the system. It turns out that a finite Zeeman splitting, A — e± — ej, changes the 
transport characteristics substantially. Figure l2"0k illustrates the gate and bias voltage 
dependence of tunnel magnetoresistance for A = 4r. First of all, a finite Zeeman 
splitting affects the TMR mainly in the regimes A (A') and B, where the second-order 
processes dominate. Furthermore, TMR exhibits a distinctively different behavior in 
the regions A and A'. It is enhanced in the case of empty dot, and reduced in the 
case of doubly occupied dot (it may even become negative in the latter case). This 
effect depends on the orientation of applied magnetic field; when the field is applied in 
the opposite direction, A = — 4r, there is an enhancement of TMR in regime A' and 
suppression in regime A, as illustrated in Fig. 120b . Thus, by changing the sign of the 
Zeeman splitting (magnetic field orientation) or applying a gate voltage, which changes 
the dot occupation, it is possible to reduce or enhance the TMR effect considerably. 
Another interesting feature of TMR displayed in Fig. I2"0k appears in the region B. 
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There is a strong asymmetry of TMR with respect to the bias reversal. The TMR is 
decreased for positive and increased for negative bias voltage. This can also be seen 
in Fig. I2"0b . but for positive transport voltage TMR is then larger than for negative 
bias voltage. The crossover between those two values of TMR takes place roughly at 
the zero bias. As a consequence, by changing the bias voltage in the small range one 
can substantially tune magnitude of the TMR effect. 

4-2. Transport through quantum dots in the Coulomb blockade regime: collinear 
magnetizations 

As we know from previous subsection, the most interesting features due to 
ferromagnetism of the electrodes, e.g. the zero-bias anomaly, occur in the blockade 
regime, where sequential tunnelling processes play a minor role. This anomaly appears 
only in the case when the dot is singly occupied at equilibrium (e CT < 0, e a + U > 0) 
and the system is in a deep Coulomb blockade regime; T, k^T <C |£ ff |, s a + U . In such 
a case the sequential tunnelling is exponentially suppressed, and cotunneling gives 
the dominant contribution to electric current [76l 11551 H"57l I171[ [T72l 1173] . Therefore, 
the sequential tunnelling can be completely ignored. In this subsection we look more 
carefully at these features taking into account only second-order contribution. 

In order to calculate the cotunneling current in the deep blockade regime one 
can employ a simplified second-order perturbation theory and master equation for the 
occupation probabilities. The rate of a cotunneling processes from lead r to lead r' , 
which change the dot state from \\) to \\'), can be written as 



with e% and st denoting the energies of initial and final states, |$*) being the state of 
the system with an electron in the lead r and the dot in state \x), whereas is a 
virtual state with e v denoting the corresponding energy. Among different cotunneling 
processes one can distinguish the single-barrier (r = r') and double-barrier (r 7^ r') 
cotunneling as well as spin-flip (\ 7^ x') an d non-spin-flip (% = \') cotunneling. The 
spin-flip processes change the spin state of the dot, whereas the non-spin-flip processes 
do not change the dot state. The current flows through the system due to double- 
barrier cotunneling processes. On the other hand, the single-barrier processes do not 
contribute directly to electric current, however, they can change the dot occupations, 
and this way also influence the current. 

The cotunneling current flowing through the system from the left to right lead is 
given by 



xx' 

where P x denotes the corresponding occupation probability. The probabilities P x can 
be found from the master equation, 




2 



(41) 




(42) 




(43) 



rr' x' 



together with the normalization condition P x = 1- 
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4-2.1. Zero-bias anomaly and its physical mechanism In order to elucidate and 
understand the anomalous behavior of TMR in the Coulomb blockade regime, we 
show in Fig. I2"la the differential conductance in the small bias regime for both parallel 
and antiparallel configurations and for several temperatures. Figure [2Tb displays the 
corresponding TMR effect. First of all, the TMR effect in the regime B for \eV\ ^> k^T 
increases with lowering temperature and approaches the Julliere's value, whereas the 
minimum at zero bias does not depend on temperature. The differential conductance 
in the parallel alignment has characteristics typical of the cotunneling regime, with a 
smooth parabolic dependence on the bias voltage. For antiparallel configuration, on 
the other hand, differential conductance has a local maximum at zero bias, followed 
by local minimum with increasing bias, as illustrated in Fig. 121b . This zero-bias 
anomaly stems from the interplay of the spin-flip and non-spin-flip single-barrier and 
double-barrier cotunneling processes [155j . The minimum in the TMR effect is a direct 
consequence of this anomalous behavior of differential conductance in the antiparallel 
configuration. 




eV/r 



Figure 21. The differential conductance (a) for the parallel and antiparallel 
configurations and the tunnel magnetoresistance (b) as a function of the bias 
voltage for different values of temperature. The maximum in differential 
conductance for antiparallel configuration at zero bias is clearly demonstrated. 
The other parameters are the same as in Fig. 1191 Figure was generated using the 
scheme for the perturbation expansion in the Coulomb blockade regime. (After 
Ref. [156] ) 

The zero-bias anomaly in the cotunneling regime is qualitatively similar to 
the anomaly due to the Kondo effect, which occurs in the strong coupling limit 
[1741 11751 1176j . There are, however, some distinct differences. First of all, processes 
responsible for the zero-bias anomaly in the cotunneling regime are of the second order 
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in tunnelling processes, while these leading to the Kondo effect are of higher order. 
The conductance in the cotunneling regime is much smaller than in the Kondo regime, 
where almost perfect transmission (G — e 2 /h) through the dot is possible owing to 
the Kondo peak in the density of states at the Fermi level. Furthermore, the Kondo 
peak occurs at temperatures lower than the so-called Kondo temperature, T < Tk, 
and exists also in the parallel configuration [5Tj 139J. 

To understand the mechanism of the zero-bias anomaly it is crucial to distinguish 
between different types of cotunneling processes: the single-barrier cotunneling 
processes shown in Fig. I2"2"b and double-barrier cotunneling processes illustrated 
in Fig. [22b . Both single-barrier and double-barrier processes can be either spin- 
flip or non-spin- flip ones. The current flows due to double-barrier cotunneling, 
whereas the single-barrier cotunneling can influence the current in an indirect way, by 
changing the spin state of the dot. In the antiparallcl configuration, there is a finite 

Spin ^ccVUPnlafinn nT1 tnA Hnt nrpcpnf cr\ in TTicr |99| ^ Tn^i rlifffront OCCUpatiOn 



a 




Figure 22. Single-barrier (a) and double- barrier (b) cotunneling processes, and 
the occupation probabilities for spin-up and spin-down electrons in the antiparallel 
configuration (c). The parameters are k^T = 0.5r, e = —U/2, U = 30r, and 
p = 0.5. (After Ref. [155]) 

probabilities for spin-up and spin-down electrons appear due to spin asymmetry 
in tunnelling processes. In equilibrium, both rates are equal and there is no spin 
accumulation, Pf = Pj. When a bias voltage is applied, and the system is in 
the antiparallel configuration, the relative amount of single-barrier cotunneling is 
diminished as compared to the double-barrier cotunneling. This is because the rate of 
single-barrier cotunneling is proportional to thermal energy, whereas that of double- 
barrier cotunneling is proportional to the bias voltage. In the nonlinear response 
regime, the magnetic state of the dot is mainly determined by the spin-flip processes 
that transfer an electron from the left to the right leads. The one shown in Fig. [2^b 
changes the dot spin from | f) to | {). Because the rate of this process is proportional 
to a product of densities of states for majority electrons, the corresponding rate is 
larger than that of the other process that changes the dot spin from | \) to | t), 
where only the minority spins are involved. This results in a nonequilibrium spin 
accumulation, Pj > P|, that increases with increasing voltage, as shown in Fig. I2"2"b. 
The initial state for the dominant spin-flip cotunneling process that contributes to 
current is | |), as sketched in Fig. I2"2"b. Thus, the conductance is diminished by spin 
accumulation. This is the mechanism by which spin accumulation gives rise to nonzero 
tunnel magnetoresistance effect, G p > G AP , see Fig. [2Tb. 
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Since the spin accumulation reduces electronic transport, any spin- flip process 
that reduces the spin accumulation should enhance the conductance. In particular, 
single-barrier spin-flip cotunneling is an example of such a process. As pointed 
above, the rate of single-barrier processes scales with U-qT while that of double- 
barrier cotunneling is proportional to max{|eV|, &bT}, which explains the mechanism 
of the zero-bias anomaly. At low bias voltage, \eV\ < fceT, single-barrier spin-flip 
processes play a significant role - they decrease the spin accumulation opening this 
way the system for the fastest double-barrier cotunneling. As a consequence, the 
current increases relatively fast with applied bias, leading to a maximum in differential 
conductance. For \eV\ ^> U-qT^ on the other hand, the relative role of single-barrier 
processes is negligible as compared to double-barrier cotunneling, and the conductance 
is reduced. Thus, the interplay between the rates of double-barrier and single-barrier 
cotunneling processes leads to the maximum in the differential conductance at the 
zero bias. 



4-2.2. Effects of spin relaxation in the dot on cotunneling transport Intrinsic spin 
relaxation in the dot can result, for instance, from spin-orbit interaction, coupling of 
the electron spin to nuclear spins, etc. However, we will not consider a particular 
microscopic mechanism of the intrinsic spin- flip processes, but simply assume that 
the spin- flip relaxation is described by a spin- relaxation time r 8 f, and is taken into 
account via a relaxation term in the appropriate master equation for the occupation 
probabilities [157] . 

0= £ {-yZ7 S P° + lZ7"Pe) - - ° % E T C , (44) 

u,i/' = L,R 

where [3 = l/ffcgT), P a denotes the probability that the dot is occupied by a spin-a 
electron, and 7^J7 CT is the cotunneling rate from lead v to lead v' with a change of the 
dot spin from a to a (a = —cr). The last term describes the spin relaxation processes, 
which in the case of spin-degenerate dot level reduces to —{Pa — P s )/t s i. 

As before, (see section 2) one can distinguish between the fast and slow spin 
relaxation limits. The former (latter) limit corresponds to the situation when the 
time between successive cotunneling events, T cot , is significantly longer (shorter) than 
the intrinsic spin relaxation time r s f . A typical spin relaxation time for quantum dots 
can be relatively long, up to ^s [1771 1178] . On the other hand, the time between 
successive cotunneling events can be estimated taking into account the fact that the 
rate of spin-flip cotunneling is generally larger than that of non-spin-flip cotunneling 
(for a finite parameter U). Assuming e-f = e[ = s, one then finds 

he 2 {e + Uf 

Tcot ~ AU 2 T 2 — ' [ ' 

with A = max{|ey|, feeT} and h = 2ttH. Assuming typical parameters [165], one can 
roughly estimate r cot to range from I0 -3 ns to I ns. 

The anomalous behavior of differential conductance in the antiparallel 
configuration results from spin asymmetry of the tunnelling processes, which gives 
rise to spin accumulation, Pf ^= P[. The intrinsic spin-flip processes in the dot reduce 
the spin accumulation and this way suppress the anomaly. This behavior is displayed 
in Fig. I23f a). where different curves correspond to different values of the parameter 
r defined as r — h/(r s fT). Thus, r = describes the case with no intrinsic spin 
relaxation, whereas the curves corresponding to nonzero r describe the influence of 
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Figure 23. The differential conductance in the antiparallel configurations (a) and 
tunnel magnetoresistance (b) as a function of the bias voltage for different spin 
relaxation r = h/(r s {T). The parameters are: k^T = 0.2F, e = — 15r, U = 30r, 
and pl — PR — 0.5. (After Ref. [157] ') 



intrinsic relaxation processes. First of all, one can note that small amount of intrinsic 
spin-flip processes enhances the zero-bias anomaly [see the curve for r = 10~ 2 in 
Fig. [23Ta)]. This is because such processes play then a role similar to that of single- 
barrier spin-flip cotunncling. In the case of fast spin relaxation, on the other hand, 
the spin accumulation is suppressed and the anomaly disappears, as can be seen in 
Fig- 1221a) for r = 1. It is also worth noting that spin-flip processes in the dot enhance 
the overall conductance in the antiparallel configuration. In the parallel configuration, 
however, the differential conductance does not depend on intrinsic relaxation. 

Modifying the conductance in the antiparallel configuration, spin relaxation 
processes suppress the tunnel magnetoresistance, as shown in Fig. 1231( b) . Since 
the zero-bias maximum in conductance is suppressed in the fast spin relaxation 
limit, the corresponding dip in TMR at small voltages is also suppressed by the 
intrinsic relaxation processes. More specifically, the dip in TMR broadens with 
increasing r and disappears in the limit of fast relaxation (see the curve for r = 1). 
An interesting feature of TMR in the presence of spin-flip scattering in the dot is 
the crossover from positive to negative values when r increases, as illustrated in 
Fig. 1123Tb). Thus, the difference between conductances in the parallel and antiparallel 
magnetic configurations persists even for fast spin relaxation in the dot, contrary 
to the sequential tunnelling regime, where such a difference disappears [147] . This 
seemingly counterintuitive behavior can be understood by taking into account the 
following two facts: (i) absence of spin accumulation in the dot for fast spin relaxation 
(Pf = Pi), and (ii) difference in the fastest cotunncling processes contributing to the 
current in the two magnetic configurations. The fastest double-barrier cotunncling 
processes involve only the majority-spin electrons of the two leads - thus, in the 
parallel configuration the fastest cotunneling processes are the non-spin-flip ones. 
They take place either via the empty-dot virtual state (for one orientation of the 
dot spin) or via the doubly occupied dot virtual state (for the second orientation 
of the dot spin). The dominant contribution to the current is then proportional to 
1/e 2 + l/(e + U) 2 . On the other hand, in the antiparallel magnetic configuration 
the fastest cotunncling processes are the spin-flip ones, which can occur only for 
one particular orientation of the dot spin. However, for this spin orientation 
cotunneling can take place via both empty and doubly occupied dot virtual states. 
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The corresponding dominant contribution to electric current is then proportional to 
[l/e-l/(e + U)} 2 = l/e 2 + l/(e + U) 2 -2/[s(e + U)}. It is thus clear that the difference 
in currents flowing through the system in the antiparallel and parallel configurations 
is equal to — 2/[e(e + U)], which results from the interference term. Since e < and 
e + U > 0, this interference contribution is positive. As a result, the current in the 
antiparallel configuration is larger than the current in the parallel configuration. 

The minimum in TMR at zero bias in the case of a symmetric Anderson model 
can be expressed as 

2p 2 (4fc B Tr 2 - e 2 h/T s{ ) 
™ Rmin " 12(1 -p 2 )k B TT 2 + (3 +p 2 )e 2 h/T s{ > (46) 
whereas for \eV\ ^> k B T and r«l one finds 

TMR - ^ [2(3-p 2 )fc B TrWVr sf ] 

* ma *~ 2(1- p 2 )(3-p 2 )k B TT 2 + (3+ p 2 )e 2 h/T sf - 

The latter formula approximates the value of TMR corresponding to the bias voltage 
at which the differential conductance has a local minimum. In the slow spin relaxation 
limit one finds TMR min = 2p 2 /(3 - 3p 2 ), and TMR max = 2p 2 /(l -p 2 ). However, in 
the limit of fast spin relaxation TMR becomes negative and is given by 

TMR min = TMR max = -2p 2 /(3 + p 2 ). (48) 



4-2.3. Effects of external magnetic field on cotunneling transport The discussion up 
to now was limited to the case of degenerate dot level. The situation changes when 
Ef =/= si, e.g., due to an external magnetic field. The level splitting is described 
by the parameter A = ej, — e-f, where the magnetic field is assumed to be along the 
magnetic moment of the left electrode. In Fig. [24] we show the bias voltage dependence 
of the differential conductance in the parallel and antiparallel configurations for 
different values of parameter r. In the limit of no intrinsic spin relaxation in the 
dot (solid line in Fig. |24|) and at low bias voltage, the dot is occupied by a spin-up 
electron and the current flows mainly due to non-spin-flip cotunneling. The spin- 
flip cotunneling processes are suppressed for |A| > |eV|,&BT, which results in the 
steps in differential conductance at |A| ~ \eV\. The suppression of spin-flip inelastic 
cotunneling was recently used as a tool to determine the spectroscopic ^-factor [165j . 
When \eV\ becomes larger than |A|, spin-flip cotunneling is allowed, consequently the 
conductance increases. However, there is a large asymmetry of differential conductance 
in the antiparallel configuration with respect to the bias reversal. To understand this 
asymmetry, it is crucial to realize that when the splitting A = e± — e-\ is larger than 
k B T, the single-barrier spin-flip cotunneling processes can occur only when the dot 
is occupied by a spin-down electron. Thus, the single-barrier processes can assist the 
fastest double-barrier cotunneling processes, but only for positive bias. This is because 
the fastest processes can occur when the dot is occupied by a spin-down electron for 
negative bias and by a spin-up electron for positive bias, leading to larger conductance 
for positive than for negative bias voltage. No such asymmetry occurs in the parallel 
configuration [see Fig. [24Tal]. as now the system is fully symmetric with respect to 
bias reversal. 

The situation changes when intrinsic spin-flip relaxation processes occur in the 
dot. For the parameters assumed in Fig.[24j the spin relaxation in the dot affects the 
conductance only for \eV\ > |A|, while for |eV| < |A| the conductance is basically 
independent of r (see Fig. [24}. This is because for |eV| < |A| and |A| » k B T, the dot 
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Figure 24. The differential conductance in the nonlinear response regime for 
different spin relaxation in the parallel (a) and antiparallel (b) configurations. 
The parameters are: k B T = 0.2T, e T = -16T, = -14T, U = 30T and p = 0.5. 
(After Ref. 1X571 1 



is predominantly occupied by a spin-up electron and the transitions to the spin-down 
state due to relaxation processes are energetically forbidden. As a consequence, the 
current flows mainly due to non-spin-flip cotunneling, irrespective of spin relaxation 
time. This scenario holds for both magnetic configurations of the system. 

When \eV\ > |A|, the spin-flip cotunneling processes can take place and the 
dot can be either in the spin-up or spin-down state. In the parallel configuration, 
Fig- [2Ua)j the conductance is slightly reduced by the spin-flip relaxation processes. 
This can be understood by realizing the fact that the fastest non-spin-flip cotunneling 
processes in the parallel configuration are more probable when the dot is occupied by 
a spin-down electron than by a spin-up one [due to smaller energy denominator, see 
Eq. ([41]) ]. Since the spin-down state (as that of larger energy) relaxes relatively fast 
to the spin- up state (which has definitely smaller energy) , this leads to a reduction in 
the conductance. On the other hand, the differential conductance in the antiparallel 
configuration is enhanced by the relaxation processes for positive bias and diminished 
for negative bias voltages. Consider first the situation for positive bias. As already 
discussed above for r = 0, an important role in that transport regime is played by 
the single-barrier spin-flip cotunneling processes, which open the system for the fast 
double-barrier cotunneling by reversing spin of the dot from the spin-down to the 
spin-up state. The relaxation processes play a role similar to that of the single-barrier 
cotunneling, and lead to a certain increase in the conductance. For negative bias 
voltage, in turn, the fast double-barrier cotunneling processes occur when the dot is 
occupied by a spin-down electron. The probability of such events is decreased by 
spin relaxation, leading to a reduced conductance. An interesting consequence of the 
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Figure 25. The differential conductance in the nonlinear response regime for 
the asymmetric Anderson model, e = —(7/4, for different spin relaxation. The 
parameters are: k B T = 0.2T, U = 60I\ p L = 0.95 and p R = 0. (After Ref. [157] 1 



enhancement (reduction) of the differential conductance for positive (negative) bias 
voltage is an increase of the asymmetry with respect to the bias reversal - see the 
curve for r = 1 in Fig. [247 b). 

4-2.4- Asymmetric situations An interesting situation occurs when the quantum dot 
is coupled asymmetrically to the left and right leads (pt, ^ pr). The differential 
conductance for a system with one electrode nonmagnetic and the other one made 
of a ferromagnet with large spin polarization (in the following referred to as strong 
ferromagnet) is shown in Fig. [55] for the case when the dot is described by an 
asymmetric Anderson model (|e| ^ e + U). 

Consider first the situation in the absence of intrinsic spin relaxation in the dot 
(solid curves in Fig. [25]) . When \eV\ ^> k^T, the influence of single-barrier cotunneling 
can be neglected |155j . The cotunneling processes which transfer charge from one lead 
to another take place via two possible virtual states - empty dot (an electron residing 
in the dot tunnels to one of the leads and another electron from the second lead 
enters the dot) and doubly occupied dot (an electron of spin opposite to that in the 
dot enters the dot and then one of the two electrons leaves the dot). Consider first 
positive bias (eV > 0, electrons flow from right to left, i.e., from normal metal to 
strong ferromagnet), and assume for clarity of discussion that the strong ferromagnet 
is a half- metallic one with full spin polarization (only spin-up electrons can then tunnel 
to the left lead). When a spin-down electron enters the dot, it has no possibility to 
leave the dot for a long time. The allowed cotunneling processes occur then via doubly 
occupied dot virtual states. In the absence of intrinsic spin relaxation in the dot, the 
only processes which can reverse the dot spin are the single-barrier cotunneling ones, 
which however play a minor role when fceT l e ^l- Thus, the current flows due to 
non-spin-flip cotunneling via doubly-occupied dot virtual states, whereas cotunneling 
through empty-dot virtual states is suppressed. The situation is changed for negative 
bias (electrons flow from strong ferromagnet to normal metal). Now, the dot is mostly 
occupied by a spin-up electron, which suppresses cotunneling via doubly-occupied dot 
virtual state and the only contribution comes from cotunneling via empty-dot virtual 
state. The ratio of cotunneling rates through the empty dot and doubly occupied dot 
virtual states is approximately equal to £ = [e/(e + U)]~ 2 . In the situation presented 
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in Fig. [2U one finds £ 3> 1. Accordingly, the conductance for negative bias is much 
larger than for positive bias voltage. 

When \eV\ becomes of the order of kpT or smaller, the rate of single-barrier 
cotunneling is of the order of the rate of double-barrier cotunneling. Therefore, the 
single-barrier processes can play an important role in transport. More precisely, single- 
barrier cotunneling processes can reverse spin of an electron in the dot and thus can 
open the system for the fast cotunneling processes. 

Intrinsic spin-flip processes in the dot have similar influence on electronic 
transport as in the case discussed in the previous section. As before, relaxation 
processes remove the asymmetry with respect to the bias reversal and suppress the 
zero-bias anomaly. Thus, the diode-like behavior can appear only in the limit of slow 
spin relaxation, and is suppressed in the limit of fast spin relaxation, as shown in 
Fig. [25] by the curves corresponding to r = 1. 



4-3. Systems with noncollinear magnetizations 

New features of transport characteristics occur when the leads' magnetizations form 
an arbitrary noncollinear magnetic configuration [1371 11401 11411 11481 [149, 150, 11611 
\TM [TSUI [TSTl [1521 [1551 \TM\ . Such a configuration can be controlled by a weak 
external magnetic field, weak enough to neglect the corresponding Zeeman splitting 
of the dot level and weaker than the effective exchange field exerted on the dot 
by ferromagnetic leads. In particular, the differential conductance is significantly 
modified by the exchange field in noncollinear configurations. This also applies to 
the TMR effect, which for arbitrarily aligned leads' magnetizations can be defined as, 
TMR = [Ip — I (<p)] 1 1 ((f) , where ip is an angle between the leads' magnetic moments. 

The system, whose spin moments Sl and Sr of the left and right lead, 
respectively form an ar' 1 ' ;. [2H The corresponding 



e+U ^ ja 




Figure 26. (color online) Schematic of a quantum dot coupled to ferromagnetic 
leads with non-collinearly aligned magnetizations. The net spin moments of the 
left Sl and right Sr lead form an angle ip. There is a symmetric bias voltage 
applied to the system. 



Hamiltonian, written for the quantization axis of the dot equivalent to that of the 
left lead, is given by Eq. ([2T>[) . where the tunnel Hamiltonian can be decomposed into 
two terms Hi — -Htl + -^tr- The term ifr/L describes tunnelling processes between 
the left electrode and the dot and takes the same form as for collinear configurations, 
while i?TR describes tunnelling processes between the dot and right lead, and acquires 
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the form 
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(49) 



The real-time diagrammatic method described at the beginning of this section 
for collinear magnetic configurations can be extended to include the effects due to 
noncollinearity of the leads' magnetic moments. The density matrix of the quantum 
dot for an arbitrary magnetic configuration is given by 

/ P ° \ 





PD = 



V 



pi pi 



T ± 1 




(50) 



/ 



The diagonal elements of the density matrix correspond to the respective occupation 

D j and Pj 



probabilities, while the off-diagonal elements P} and Pj describe the dot spin S, with 



S x = ReP| , = TxaPj, and 5 Z = ^P| — P£J /2. The density matrix elements can 

be determined from the corresponding kinetic equation, which in the steady state and 
for spin-degenerate dot level can be written as [1041 11611 11661 1167] 







EpXi^XiXl 
X' 2 X2X2 

x'i,x' 3 



(51) 



By expanding the self-energies and density matrix elements, one can calculate the 
sequential and cotunneling current. 

In the following we focus on transport in the Coulomb blockade regime, |e|, |e + 
U\ ^> r, ftsT, and when the dot is singly occupied, e < < e+U. ^From the discussion 
above we know that there is an anomalous behavior of the differential conductance in 
the small bias regime for antiparallel magnetic configuration. Now, we will consider 
how this anomaly changes when magnetic configuration varies continuously from 
antiparallel to parallel alignment. 



4-3.1. Symmetric Anderson model We consider first the symmetric Anderson model, 
e = —U/2. The exchange field vanishes then for an arbitrary magnetic configuration, 
and both differential conductance and TMR change monotonically when going from 
antiparallel to parallel magnetic configurations. However, for e ^ —U/2, the effects 
of exchange field become important and lead to nontrivial behavior of transport 
characteristics, as will be discussed later. 

The differential conductance G as a function of the bias voltage is shown in 
Fig- l27T a) for several values of the angle <p. The conductance varies monotonically 
with the angle between leads' magnetic moments. When the angle increases from zero 
to 7r, the anomaly emerges at small values of (p and its relative height increases with 
increasing angle, reaching a maximum value at (p = tt. 

The exchange field for a symmetric Anderson model is negligible, and the average 
dot spin tends to zero in the linear response regime. The angular dependence of the 
linear conductance can be then expressed as 



G 



2he 2 



3-P 2 1 



2 sin 2 - 



(52) 
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Figure 27. (color online) Differential conductance G = dl/dV (a) and tunnel 
magnetoresistance TMR (b) as a function of the bias voltage V for indicated 
values of the angle ip between magnetic moments, calculated for symmetric 
Anderson model. The parameters are: k-^T = 0.5r, e = — 15F, U = 30r, and 
p L = p R EE p = 0.5. (After Ref. [TSU I 



The variation of G with <p at low bias is thus characterized by the factor 1 + 2 sin 2 (</3/2), 
which leads to maximum (minimum) conductance in the parallel (antiparallel) 
magnetic configuration. Such behavior is typical of a normal spin- valve effect. 

The bias dependence of the associated TMR is shown in Fig. l2W b) for several 
values of the angle (p. The zero-bias anomaly in the differential conductance [see 
Fig. l2"TTa)] leads to the corresponding anomaly (dip) in TMR at small bias voltages. 
The dip in TMR decreases when magnetic configuration departs from the antiparallel 
alignment, and eventually disappears in the parallel configuration. The variation 
of TMR with the angle is monotonic, similarly as the angular variation of the 
differential conductance. The dependence of TMR on the angle ip at zero bias is 
given approximately by the formula 

2p 2 sin 2 f , x 

TMR = 2 „ „. . 53) 

3-p 2 (l + 2sin 2 f) V ' 

Now, the angular dependence of TMR is governed by sin 2 (<y2/2), which gives maximum 
TMR in the antiparallel configuration and zero TMR in the parallel one. 



4-3.2. Asymmetric Anderson model The transport characteristics, described above 
for the symmetric Anderson model, are significantly modified when the model becomes 
asymmetric, i.e. when e ^ —U/2. This can be realized by shifting the dot level position 
by a gate voltage applied to the dot. Different contributions to the exchange field do 
not cancel then, and the resulting effective exchange field becomes nonzero. And 
this exchange field has a significant influence on transport properties. The strength 
of effective exchange field grows with deviation of the Anderson model from the 
symmetric one, described quantitatively by 2e + U (with 2e + U = for the symmetric 
model). 

The leading contribution to the exchange field comes from the first-order diagrams 
of the perturbation expansion in terms of the real-time diagrammatic technique, 
whereas the current flows due to the second-order tunnelling processes. This leads 
to two different time scales which determine transport characteristics. Transport 
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Figure 28. (color online) The differential conductance G = dl/dV (a) and tunnel 
magnetoresistance TMR (b) as a function of the bias voltage V for indicated values 
of the angle ip for asymmetric Anderson model e = — 12r and U = 30I\ The other 
parameters are the same as in Fig, 1271 (After Ref. I161p 



properties are thus a result of the interplay of the first- and second-order processes. 
Although the first-order processes do not contribute directly to current, they influence 
transport via modification of the dot spin. 

The differential conductance is shown in Fig. I^ST a) as a function of bias voltage 
for an asymmetric Anderson model, e ^ —U/2, and for different values of (p. The 
low-bias differential conductance becomes enhanced in a certain range of the angle ip, 
leading to a nonmonotonic dependence of the conductance G on the angle between 
the leads' magnetizations. This nonmonotonic behavior is also visible in TMR, see 
Fig.HSSTb). There is also a range of the angle ip, where TMR changes sign and becomes 
negative in the small bias region, i.e. the corresponding conductance is larger than 
that in the parallel configuration. 

The most characteristic features of transport characteristics in the presence of 
exchange field are the enhanced differential conductance at a non-collinear alignment 
and its rapid drop when the system approaches the antiparallel configuration. The key 
role in this behavior is played by the first-order processes giving rise to the exchange 
field. These processes lead to the precession of spin in the dot, which facilitates 
tunnelling processes and leads to an increase in the conductance as compared to that 
in the parallel configuration. When the configuration becomes close to the antiparallel 
one, the first-order processes become suppressed and the conductance drops to that 
for antiparallel alignment. The nonmonotonic behavior of the differential conductance 
with cp leads to a nonmonotonic dependence of TMR. The effects due to exchange 
field give rise to a local minimum in TMR at a non-collinear magnetic configuration. 
Moreover, in this transport regime TMR changes sign and becomes negative. When 
the magnetic configuration is close to the antiparallel one, TMR starts to increase 
rapidly reaching maximum for ip = n. The negative TMR and its sudden increase 
when the configuration tends to the antiparallel one are a consequence of the processes 
leading to nonmonotonic behavior of the differential conductance, as described above. 

To understand more intuitively the above presented behavior of the differential 
conductance and TMR at low bias voltage and close to the antiparallel configuration, 
one should consider two different time scales. One time scale, r proc , is established by 
the virtual first-order processes responsible for the spin precession due to exchange 
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field, 

(54) 



IVec| ~ -^rP sin ip In 



U 



The second time scale, r cot , is associated with second-order processes which drive the 
current through the system. At low temperature and low bias voltage, the cotunneling 
rate can be expressed as 

p2 \eV\U 2 
7-cot ~ ih ( - 1+P ^ 1 ' pCOSlf K 2 (e + U) 2 ' (55) 

The rate r co t depends linearly on the applied voltage, whereas r prcc is rather 
independent of V. As a consequence, at low bias and for non-collinear configuration, 
the exchange field plays an important role leading to a nonmonotonic dependence 
of differential conductance on the angle between the leads' magnetizations. When 
magnetic configuration is close to the antiparallel one, the spin precession rate is 
deceased (|r^ c | ~ siny) and, at certain angle, the rate of spin precession becomes 
comparable to the cotunneling rate. This gives rise to a sudden drop (increase) 
in differential conductance (TMR). We note that the nonmonotonic dependence of 
differential conductance and magnetoresistance has also been observed in quantum 
dots in the strong coupling limit [1401 1141] . 



5. Transport through multi-level quantum dots connected to 
ferromagnetic leads 

Up to now we discussed theoretical aspects of electronic transport through the simplest 
quantum dots, i.e. the dots with only one orbital level. Real dots are usually more 
complex and their electronic spectrum includes many orbital levels taking part in 
electronic transport. This may significantly change transport characteristics and also 
may lead to qualitatively new features pIU [1861 QUI HM HH [190]. In this section 
we will consider some of the new effects in transport through multi-dot systems. 



5.1. Sequential transport in two-level dots 

The schematic of a two-level quantum dot coupled to ferromagnetic leads with collinear 
magnetizations is shown in Fig. 1291 The quantum dot is described by the following 
Hamiltonian 

H D = ^2 £ i n i° + U Y1 n ^ n n +U 'Y1 nxvtiw - — ^ (n iT - n jl ) , (56) 

3" 3 "a 1 j 

where n JCT is the particle number operator, nj a — dj^dja, d^ a (dj a ) is the creation 
(annihilation) operator of an electron with spin a on the jth level (j = 1,2), and e 3 - 
is the corresponding single-particle energy. The on-level Coulomb repulsion between 
two electrons of opposite spins is described by U , whereas the inter-level repulsion 
energy is denoted by U'. The forth term in Eq. ([56")) describes the Zeeman energy, 
with A = gfJ-sB being the Zeeman splitting of the energy levels (B is an external 
magnetic field along the magnetic moment of the left electrode). To present the main 
features of transport characteristics, it is convenient to introduce the level spacing 
Se = E2 — E\ and also define E\ = e. 

In the following we analyze the current /, differential conductance G, and the 
Fano factor F in the parallel and antiparallel magnetic configurations, as well as 
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Figure 29. (color online) Energy diagram of a two-level quantum dot coupled 
to ferromagnetic leads. For clarity reasons the energy diagram is shown here for 
A = U' = 0. The leads' magnetizations can form either parallel or antiparallel 
configurations. The arrows indicate the net spin of the leads. 



the corresponding TMR for two-level quantum dots. The Fano factor, F — S/S p , 
describes the deviation of the zero-frequency shot noise S from the Poissonian shot 
noise S„ — 2e\I\. The nresented results have been obtained within the real-time 
diag, □ □ □ 
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Figure 30. (color online) The differential conductance G = dl/dV as a function 
of the bias voltage and level position in the parallel magnetic configuration (a) 
and TMR (b) for the parameters: k B T = T, 5s = 25I\ U = 501", A = 0, 
p L = p R = p = 0.7, and T rj = T/2 (r = L,R, j = 1,2). (After Ref. (155) 1 

5.1.1. Quantum dots symmetrically coupled to ferromagnetic leads Typical variation 
of the differential conductance with the bias voltage V and level position s (gate 
voltage) is shown in Fig. l5DT a) for the parallel magnetic configuration. Conductance 
in the antiparallel configuration is qualitatively similar to that in the parallel 
configuration, but generally smaller, which leads to nonzero TMR effect. The 
diamonds in Fig. (30T a) around V = correspond to the Coulomb blockade regions. 
When lowering position of the dot levels, the charge of the dot changes successively. 
The dot is empty for e > 0, occupied by one electron for > e > —U, doubly occupied 
for —U > e > — (2U+5e), occupied by three electrons for — (2U +5e) > e > ~{2>U +5e), 
and the two orbital levels of the dot are fully occupied for —(3U + Se) > e. In all these 
transport regions the dot is in a well-defined charge state, and the sequential tunnelling 
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Figure 31. (color online) Current (a) in the units of io = eT/h and Fano factor 
(b) in the parallel (solid line) and antiparallcl (dashed line) configurations as well 
as tunnel magnetoresistance (a) as a function of the bias voltage for e = U/2. 
The other parameters are the same as in Fig. 1301 (After Ref. [1880 

is exponentially suppressed. If the bias voltage is increased above a certain threshold 
voltage, the current flows due to first-order (sequential) tunnelling processes. When 
the thermal energy is low enough, one observes then a well-resolved step in the current 
as a function of the bias voltage. In the density plots shown in Fig. [50b . this can be 
seen in the form of lines that clearly separate the Coulomb blockade regions from 
transport regions associated with consecutive charge states taking part in transport. 
When the bias voltage increases further, additional steps [and consequently lines in 
Fig. 1301 a)] arise at voltages where new states becomes active in transport. 

The corresponding TMR as a function of the bias and gate voltages is shown in 
Fig. l3"07 bV It is worth noting that TMR takes now several well-defined values. Such 
behavior of TMR is significantly different from that for a single- level quantum dot, 
where TMR in the sequential tunnelling regime acquires only two values [156] . As in 
single-level dots, TMR in the linear response regime is independent of the gate voltage 
and is given by p 2 /{l~p 2 ), i.e. half of the TMR in the Julliere model [13]. Fi gurel30l 
also shows that, when increasing the bias voltage V and keeping constant position of 
the dot levels, the current and TMR acquire some specific and well-defined values in 
different transport regions. As shown in Ref. [188j . the current and TMR at these 
plateaus can be approximated by simple analytical formulas. 

This behavior is presented in more details in Fig. I3l7 a). where the bias voltage 
dependence of the current and TMR is shown explicitly for the case when the dot 
level is above (e = U/2) the Fermi level of the leads at equilibrium. The current in 
both magnetic configurations and the associated TMR exhibit characteristic plateaus 
which correspond to different transport regions. The corresponding Fano factors Fp 
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and _Fap i n the parallel and antiparallel magnetic configurations are shown in Fig. 131b . 
Similarly as current and TMR, the Fano factor acquires roughly constant values, 
different in different transport regions. Due to the spin asymmetry in the coupling 
of the dot to external leads, the bias dependence of the Fano factor is significantly 
different from that in the corresponding nonmagnetic situations |185j . For both 
magnetic configurations of the system, the Fano factor depends on the polarization 
factor p (differently in the two configurations, in general). If \eV\ <C k^T, the Fano 
factor becomes divergent due to the thermal noise, which dominates in this transport 
regime; in the case of V = 0, the noise is given by S — Ak^TG 1 " 1 , with G lm being 
the linear conductance, leading to a divergency of the Fano factor |116[ 1191] , Apart 
from this, in some transport regions we find Fp > Fap, while in the other ones 
Fp < Fap ■ The ratio Fp /Fap depends generally on the spin polarization of the leads 
p |188j . Furthermore, we note that if the leads are half- metallic, the Fano factor in 
the parallel configuration diverges as p — > 1. This increase of the Fano factor is due to 
the enhanced spin asymmetry in transport processes through the dot |146l 11921 1193j . 
On the other hand, in the antiparallel configuration the Fano factor tends to unity 
for p — > 1, except for the Coulomb blockade regime with two electrons trapped in the 
dot, where the shot noise is super-Poissonian. 

5.1.2. Quantum dots asymmetrically coupled to the leads When one of the leads 
is half- metallic (p = 1) and the other one is nonmagnetic (p — 0), transport 
characteristics become asymmetric with respect to the bias reversal. Furthermore, the 
current can be suppressed in certain bias regions, and this suppression is accompanied 
by the occurrence of NDC. This basically happens when the electrons residing in the 
dot have spin opposite to that of electrons in the half-metallic drain electrode. In 
Fig. [321 we show the current and Fano factor for a quantum dot coupled to half- 
metallic (left) and nonmagnetic (right) lead as a function of the bias voltage. For 
the situation shown in Fig. 1321 current is suppressed in certain regions of positive bias 
voltage. As follows from Fig. l327 a). there are three such blockade regions, labelled with 
the consecutive numbers. On the other hand, for negative bias voltage, the current 
changes monotonically with the transport voltage, as shown in the inset of Fig. 1327 a). 

The mechanism of the blockade can be described as follows. In the blockade 
region (1), 1 < eV/U < 2, the dot is in the state | |)|0) and the spin-down electron 
residing in the dot has no possibility to tunnel further to the left lead, which leads to 
suppression of the current. Here, the first (second) ket corresponds to the first (second) 
orbital level of the dot. The blockade region (2), 2 < eV/U < 3, is associated with 
the occupation of the states | |)|0) and |0)| J.). The current is then prohibited due to 
the full occupation of the single-particle spin-down states. When increasing the bias 
voltage further, 3 < eV/U < 4, the blockade of the current becomes suppressed [see the 
plateau between regions (2) and (3) in Fig. 1327 a)]. which is due to a finite occupation 
of state | Tl)|0)- Although tunnelling of spin-down electrons is then blocked, the 
current is still carried by spin-up electrons. In turn, the blockade region (3) occurs 
for 4 < eV/U < 5, where the dot is in the triplet state | |)| J.) and tunnelling is 
also suppressed. Thus, the current is blocked when the total dot spin S z is either 
S z = —1/2 or S z = —1, i.e., the spin of electrons on the dot is opposite to that of 
electrons in the half-metallic lead. There is no suppression of current for negative 
voltage. This is because the electrons residing in the dot can always tunnel to the 
nonmagnetic drain electrode. It is also interesting to note that a pure triplet state is 
formed in the region (3) [US [194]. 
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Figure 32. The current (a) in units of Iq = eT/h and Fano factor (b) as a 
function of the bias voltage for pl = V = 0.95, 0.98, pr, = 0, while the other 
parameters are the same as in Fig. I30I The inset in part (a) shows the current in 
the whole range of the bias voltage. (After Ref. |188p 

The current blockade in the regions (1) to (3) in Fig. I3"2f a) for positive bias is 
not due to the charging effects as in the Coulomb blockade regime, but due to a 
particular occupation of the dot spin state. Such blockade is frequently referred to 
as the Pauli spin blockade, and has already been found in single-dot and double-dot 
systems [TM1 WM [T37I UM[ UM\ . It is worth noting that the current blockade in 
Fig. I3"!?l a) is not complete - there is a small leakage current in each blockade region, 
which results from the fact that the assumed spin polarization of the half-metallic lead 
is not exactly equal to unity. When spin polarization is increased, the current in the 
blockade regions decreases (compare the curves for p = 0.95 and p — 0.98). 

The spin blockade of charge current leads to the super-Poissonian shot noise, i.e., 
the corresponding Fano factor is larger than unity. On the other hand, the Fano factor 
outside the spin blockade regions is sub-Poissonian (smaller than unity), as shown in 
Fig- EUpj). The enhancement of the shot noise in the Pauli blockade regions is a 
consequence of large spin asymmetry in the tunnelling processes. The occurrence of 
a spin-down electron on the dot prevents further tunnelling processes for a longer 
time, while spin-up electrons on the dot can escape much faster, allowing further 
tunnelling processes. This bunching of fast tunnelling processes gives rise to large 
current fluctuations, and consequently also to Fano factors much larger than unity. 

Recently an interacting three-terminal quantum dot with ferromagnetic leads was 
considered in Ref. |192[ 1193] . The dot operated as a beam splitter - one contact was 
a source and the other two acted as drains. The authors found a dynamical spin 
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blockade (spin-dependent bunching of tunnelling events) and positive zero-frequency 
cross-correlations of the current in the drain electrodes. 

5.2. Cotunneling in two-level quantum dots 

In all transport characteristics discussed above only sequential tunnelling processes 
were taken into account, while higher order processes, in particular the cotunneling 
ones, have been neglected. From our discussion on single-level dots we already 
know that higher order contributions are particularly pronounced in the blockade 
regions, where they have a significant influence on TMR and Fano factor. The role 
of cotunneling processes in two-level dots was studied in a recent paper [200) . It has 
been shown there that the cotunneling processes lead to a significant enhancement 
(or reduction) of TMR in the blockade regime. Outside the blockade regions, 
however, TMR is determined mainly by sequential transport, so the contribution 
from cotunneling processes is rather minor. Similarly, the cotunneling processes 
significantly modify the Fano factor in the blockade regions, while outside these regions 
the Fano factor is only weakly sensitive to the cotunneling. More specifically, the 
Fano factor in the blockade regions is reduced in comparison to that obtained in 
the sequential transport limit. However, it is still larger than one, indicating super- 
Poissonian shot noise in the blockade regime due to bunching of the fastest cotunneling 
processes. 

There is also another interesting phenomenon in the blockade regime, which 
resembles the zero bias anomaly observed in single-level dots in the antiparallel 
magnetic configuration. This takes place when each orbital level of the dot at 
equilibrium is occupied by a single electron, and the system is in the Coulomb blockade 
regime |189j . The Hamiltonian H of the system is given by Eq. where £/d 

additionally includes the exchange term and can be expressed by the formula [20111202] 

H D = ^ SjUja + U 2J n^Uji +U'^2 n l<r n 2a' - J d^di^dl^s^afj^s . (57) 
jtr j aa' a0^f6 

The last term in Hp corresponds to the exchange energy according to the Hund's rule, 
with J being the respective exchange coupling and a denoting a vector of Pauli spin 
matrices. 

Six different two-particle states of the dot are possible, these are: three singlets 
\S = 0,M = 0>x = (| T)| I) - I 1)1 T»/V2, |0,0> 2 = I U)|0), |0,0> 3 = |0)| U>, and 
three triplets |1,0) = (| T)| I) + I 1)1 T))/A |1,1) = | T)| T) and |1,-1) = | |)| |). 
In the case of finite level spacing, Se > k^T, Y, and J < 5e, the lowest singlet state is 
|0,0)2. Transport characteristics strongly depend on the ground state of the system. 
It is therefore useful to introduce the difference between the energy of the lowest lying 
singlet (es) and triplet (ex) states, Ast = £s — £t = J — Se. 

For Ast < 0, the ground state of the dot is a singlet, |0, 0)2, whereas for Ast > 0, 
the ground state is a triplet, which is three-fold degenerate, |1, 0), |1, 1), |1, —1). On 
the other hand, for Ast = 0, the dot is in a mixed state and the occupation of singlet 
and each triplet is equal at equilibrium and given by 1/4. 

Figure[33]presents the bias voltage dependence of differential conductance G in the 
parallel and antiparallel magnetic configuration and for several values of Ast = J—Se. 
Let us consider the situation in Fig 1331c) , which corresponds to degenerate singlet 
and triplet states, Ast = 0. The differential conductance in the parallel configuration 
exhibits a smooth parabolic dependence on the bias voltage, whereas in the antiparallel 
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Figure 33. (color online) The differential conductance in the parallel (solid line) 
and antiparallel (dashed line) magnetic configuration of the system as a function of 
the bias voltage for several values of Agx = J — Se as indicated: (a) Agx = — 51", 
(b) A ST = -2.5r, (c) A ST = 0, (d) A ST = 2.51". The other parameters are: 
k B T = Q.5T, e = -60T, Se = 5V, U = V = 40r, and p = 0.5. (After Ref. 1X89) 1 



configuration there is a maximum at the zero bias voltage. This effect bears a 
resemblance to that found in the case of singly occupied one-level quantum dots 
discussed in previous sections [155 . Here, however, the mechanism leading to the 
maximum is different - the zero-bias peak appears due to cotunncling through the 
singlet and triplet states of the dot. In the case of A§t = and at low bias voltages, 
all the four dot states, i.e. |0, 0)2, 1 1 , 0} , |1,— 1), |1,1), participate in transport on 
an equal footing. Consequently, the current flows due to both spin-flip and non-spin- 
flip cotunneling processes. To understand the mechanism leading to the zero-bias 
peak, one should bear in mind that in the antiparallel configuration the spin-majority 
electrons of one lead tunnel to the spin-minority electron band of the other lead. 
For example, for positive bias voltage (electrons tunnel then from the right to left 
lead), the spin-| electrons can easily tunnel to the left lead (the spin-| electrons are 
the majority ones), while this is more difficult for the spin- J, electrons (they tunnel 
to the minority electron band). Thus, the occupation of state 1 1 , — 1) (| J.)| J.)) is 
increased with increasing bias voltage, while the occupation of state |1,1) (| T)| T)) 
decreases. The unequal occupations of these triplet states lead to a nonequilibrium 
spin accumulation in the dot, -P|i.i) — P\i,-i) < 0. It is further interesting to note 
that in the antiparallel configuration the possible non-spin-flip cotunncling processes 
are proportional to T^r^ and , whereas the spin-flip cotunneling is proportional 

to r^r^ and r^r^. It is clear that the fastest cotunneling processes are the ones 
involving only the majority spins, i.e. T^. However, because of nonequilibrium spin 
accumulation, with increasing the bias (V > 0), the dot becomes dominantly occupied 
by majority electrons of the right lead, Pn—i) > P\i,i), an d the processes proportional 
to r^F^ are suppressed. As a consequence, the differential conductance drops with 
the bias voltage, leading to the zero-bias peak. This is thus the nonequilibrium 
accumulation of spin 5* = 1 which is responsible for the maximum in G at low bias 
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voltage. Because the dot is coupled symmetrically to the leads, there is no spin 
accumulation in the parallel configuration and the differential conductance exhibits a 
smooth parabolic dependence. 

The maximum in the conductance at zero bias disappears for Ast < as well as 
for Ast > 0. Instead of the narrow maximum, a relatively broad minimum develops 
then in the conductance for both magnetic configurations, as shown in Fig l337 a.b) or 
for one magnetic configuration as in Fig l33T d). Consider for instance the case Ast < 0, 
shown in Fig. l33Ta.b). At low bias voltage the dot is occupied by two electrons on 
the lowest energy level and the ground state is singlet, S — 0. Current is mediated 
then only by non-spin-flip cotunneling processes. Once |eV| > |Ast|, the triplet states 
start participating in transport leading to an increase in the differential conductance 
at \eV\ w | Ast |- Thus, the suppression of cotunneling through 5=1 states gives rise 
to a broad minimum in the differential conductance. The width of this minimum is 
determined by the splitting between the singlet and triplet states, 2|Ast|- 
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Figure 34. (color online) Tunnel magnetoresistance as a function of the bias 
voltage for Agx = —5, —2.5, 0, 2. 51", as indicated in the figure. The other 
parameters are the same as in Fig. 1331 (After Ref. [189] ^ 



The bias voltage dependence of TMR is shown in Fig. [M] Generally, TMR 
displays a nontrivial dependence on the ground state of the dot. When the dot is 
occupied by a singlet, Ast < 0, TMR displays a maximum plateau at low bias, which 
can be approximated by, TMR S=0 = 2p 2 /(l - p 2 ) [155] . This is the Julliere value of 
TMR. Here, it results from the fact that in this transport regime current flows due 
to non-spin-flip cotunneling. On the other hand, for Ast > 0, the TMR exhibits a 
minimum at zero bias, which can be seen in Fig. [M] for Ast = and Ast = 2.5I\ 



5.3. Multi-level quantum dots based on single-wall carbon nanotubes 

When the number of discrete levels in the dot (relevant for electronic transport) 
becomes larger than two, transport characteristics of the dots attached to 
ferromagnetic leads become more complex and reveal further interesting features. Such 
features were observed experimentally mainly in molecule-based quantum dots. When 
a large natural molecule is weakly attached to metallic leads, it can be treated simply 
as a multi-level QD [203, 204, 205]. As a specific case, we discuss transport properties 
of a single-wall metallic carbon nanotube (CNT) weakly coupled to ferromagnetic 
leads ^[2Tm[2TM[27M[2T0] . 

Sequential transport in CNTs weakly coupled to nonmagnetic and ferromagnetic 
leads was considered in a recent paper |206j , where TMR and the Fano factor have been 
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Figure 35. (color online) The differential conductance in the parallel 
configuration (a) and tunnel magnetoresistance (b) as a function of bias and 
gate voltages. The parameters are: A = 8.4 meV, U/A = 0.26, J/A = 0.12, 
SU/A = 0.04, S/A = 0.1, k B T/A = 0.025, p L =p R = 0.5, x = 0.14, and T = 0.2 
meV. 



calculated within the real time diagrammatic technique. The system was modelled by 
the Hamiltonian in the form introduced by Oreg et al. |211j . 
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(58) 



where n 



d ■ 



and d\ ja (d 



j, w<j ) is the creation (annihilation) operator of an 
electron with spin a on the Jth level in the subband \x (/i — 1, 2). The corresponding 
energy of the jth discrete level in the subband n is given by e^j = j A + (/x — 1)5, 
where A is the mean level spacing and S describes the energy mismatch between the 
discrete levels corresponding to the two subbands. The second term in Eq. (|5"8|) stands 
for the electrostatic energy of a charged CNT, with U denoting the charging energy 
and No being the charge on the nanotube induced by gate voltages. The next term 
corresponds to the on-level Coulomb interaction with SU being the relevant on-site 
Coulomb parameter. Finally, the last term in Eq. (|58|) describes the exchange energy, 
with J being the relevant exchange parameter. 

In Fig. [35] we show the differential conductance in the parallel configuration and 
associated TMR, calculated in the first order approximation. The results are for a 
specific choice of parameters, namely 6U + J > 5. The latter condition indicates 
that the following sequence of the ground states is realized when filling the CNT with 
electrons while sweeping the gate voltage (for V = 0): 5 = 0, 1/2, 1, 1/2, where S is 
the total spin of the nanotube. This means that in a certain transport regions, the 
systems is in the triplet state at equilibrium. Fig. [35T a) clearly reveals the fourfold 
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pattern of the conductance spectra, associated with filling of the consecutive levels in 
the two electron subbands of the CNT. The corresponding TMR is shown in Fig.l^STb - ), 
and reveals a complex variation with both gate and bias voltages. 

As shown in Ref. 206J, shot noise calculated in the first order approximation is 
super-Poissonian in the blockade regions, and sub-Poissonian in the transport regions 
were sequential tunnelling dominates. The corresponding Fano factor in the latter 
case is close to 1/2, while in the former one is above 1. Detailed calculations taking 
into account second order (cotunneling) processes show that the Fano factor in the 
blockade regions is generally larger than 1, although it is reduced in comparison to 
that derived in the first order approximation, similarly as in the case of two-level dot. 

6. Kondo effect in quantum dots coupled to ferromagnetic leads 

When a quantum dot or a large molecule is strongly coupled to external leads, low- 
temperature transport characteristics reveal features that are typical of the Kondo 
phenomenon [ED EH QH [1751 HZSl [212] . To analyze the Kondo effect [213], one may 
use various theoretical techniques. For instance, one can use the perturbative methods, 
like the real-time diagrammatic technique discussed in previous sections. However, 
one has to go beyond the approximations used in the discussion of the sequential 
and cotunneling transport. Another approximate method which is frequently used to 
study electronic transport through quantum dots, including the Kondo regime, is the 
equation of motion method for the relevant Green function. Alternatively, one may 
employ nonperturbative methods, such as for example the numerical renormalization 
group method [2T31 [2151 I2IS1 I2T7] . 

The Kondo effect [78] in electronic transport through a quantum dot has been first 
predicted theoretically in Refs. [218, 219J. Since then, it has become well documented 
experimentally pi! [220] [2211 [2221 [223]. The dot with an odd number of electrons 
has a local spin which at low temperatures, k^T < fceTk <C T, and in the presence of 
strong coupling to the electrodes behaves effectively like a magnetic Kondo impurity. 
Exchange coupling with the leads' electrons gives rise to spin fluctuations in the dot 
and screening of the dot's spin. This, in turn, leads to the formation of a singlet 
state and Kondo resonance peak at the Fermi level in the dot's density of states. This 
Kondo-correlated state leads to an increased transmission through the dot, and also 
gives rise to a sharp zero-bias anomaly in the current-voltage characteristics. The 
successful observation of the Kondo effect in semiconductor-based quantum dots, as 
well as in molecular quantum dots based for instance on carbon nanotubes [224, 223] or 
other single molecules [226, 227 attached to metallic electrodes opened new possibility 
to study the influence of ferromagnetism on the Kondo phenomenon [44], [5H 1228] . 

In this section we briefly discuss how the Kondo effect is modified by 
ferromagnetism of the electrodes. First, we note that in the extreme case of half- 
metallic leads, i.e. in the absence of minority-spin electrons, the screening of the dot 
spin is not possible in the parallel magnetic configuration. Consequently, no Kondo- 
correlated state can be then formed in this particular magnetic configuration. 

The possibility of the Kondo effect in a quantum dot attached to ferromagnetic 
electrodes was discussed in a number of papers [1371 USSl ESI EH [2291 ESDI E2H [2321 
233, 234, 235J. It was shown, that the Kondo resonance in the parallel configuration 
is split and suppressed in the presence of ferromagnetic leads [1391 12311 12321 12351 1234) . 
However, it was also demonstrated that this splitting can be compensated by an 
appropriate external magnetic field which restores the Kondo effect [139| 1231] . 
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Let us begin with the collinear (parallel and antiparallel) magnetic configurations. 
For the antiparallel configuration and vanishing magnetic field and bias voltage, the 
model is equivalent to a quantum dot coupled to a single lead with density of states 
PL] + PR] = Phi + Pm [2181 1219j . Thus, the usual Kondo resonance forms in such a 
case, which is the same as for nonmagnetic electrodes [78]. The situation changes in 
the parallel configuration, where there is an overall asymmetry for up and down spins, 
say pt\ + PR] > pli + Pri- This significantly reduces (or even suppresses) the Kondo 
effect. 



6.1. Poor man's scaling approach 

Some basic relations concerning the Kondo effect, and particularly those concerning 
the Kondo temperature, can be derived from a simple poor man's scaling [236 
performed in two stages [237] . First one reduces the energy scale of the effective 
electron bandwidth D starting from Dq. This leads to a renormalization of the level 
position edcr according to the scaling equations [139] 

dln(D /D) ~ 2tt ' [ ' 

This leads to a spin splitting of the level, which in the presence of a magnetic field 
simply adds to the initial Zeeman splitting Ae d . As a result, one finds 

A? d - ? dT - ? di = -(l/7r)Pr ln(D /D) + Ae d . (60) 

The scaling of Eq. ([59]) is terminated at D ~ — [237 . 

The strong-coupling limit can be reached by tuning the external magnetic field B 
in such a way that the total effective Zeeman splitting vanishes, Aed = 0. In the second 
stage [236] . spin fluctuations are integrated out. The Schriefier- Wolff transformation 
[78] allows to map the Anderson model, with renormalized parameters D and e d , to 
the effective Kondo Hamiltonian [139] 

#Kondo = J+S + ^ a tkl a r'q] + a t q ] a r'kl 
rr'kq rr'kq 

+ S Z I J z] ^ a lq] a r'q'1 ~ J *l a lkl a r'k'l , (61) 

y rr'qq' rr'kk' I 

with J + = J_ = J 2 | = J z i = |T| 2 /|ed| = Jo m the large-C/ limit. Although initially 
identical, the three coupling constants J + = J_ = J±, J z |, and J z ^ are renormalized 
differently during the second stage of scaling. The relevant scaling equations are 



' /(/ ' ' ' =P±J±{P]J z ]+PiJ zi ) (62) 



d\n{D/D) 

d(pcrJza) 



dln(D/D) 21PM2 (63) 

All the coupling constants reach the stable strong-coupling fixed point J± = J z ^ = 
J z l = oo at the Kondo energy scale, D ~ ksTx- For the parallel configuration, 
the Kondo temperature in leading order depends on the polarization p in the leads 
following the formula 
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The Kondo temperature is maximal for nonmagnetic leads, p = 0, and vanishes for 
p — > 1. The unitary limit for the parallel configuration can be reached by tuning the 
magnetic field, as discussed above. In this case, the maximum conductance through 
the quantum dot is the same as for nonmagnetic leads, i.e. G^ ax a = e 2 jh per spin. 

The scaling procedure was extended in Ref. [1441 11451 1238) to noncollinear 
configurations, where instead of Eq. (|64|) one finds 

K ' I (PT+PiVo pcos(0/2) / 1 ' 

This formula shows explicitly how the Kondo temperature varies with the angle 9 
between magnetic moments of the leads. 



6.2. Numerical renormalization group 

Numerical renormalization-group (NRG) technique [78l I214[ I217j is one of the 
most powerful and accurate methods available currently to study strongly-correlated 
systems in the Kondo regime. This technique has been adapted recently to the case 
of a quantum dot coupled to ferromagnetic leads [1451 12311 1232j . A simple way to 
model the ferromagnetic leads in the standard NRG procedure is to take the density 
of states in the leads to be constant and independent of electron spin, p rcr (oj) = p, 
the bandwidths to be equal = Dt, and lump all the spin-dependence into the 
spin-dependent coupling parameter, T r<7 (tj), which can be taken as independent of 
energy. 

The NRG method, with recent improvements related to high-energy features and 
finite magnetic field [2391 GUI BUI EE], is a well-established method to study the 
Kondo impurity (quantum dot) physics. Using this technique one can calculate the 
level occupation n^a = (d^da), and the spin-dependent single-particle spectral density 
A a (uj) for arbitrary temperature T, magnetic field B and polarization p. Using the 
spectral function one can find the spin-resolved linear conductance 

r e 2 r Lff r Rg r°° 9/M 

with f(uj) denoting the Fermi function. 



6.3. Nonequilibrium transport: zero-bias anomaly 

The Kondo effect is usually observed as a zero bias anomaly in the differential 
conductance. This requires theoretical methods applicable to nonequilibrium 
situations. One of such techniques is based on the nonequilibrium Green functions. 
The latter are usually determined from the relevant equation of motion and the 
appropriate decoupling scheme. To describe the main features of the Kondo effect 
in transport through quantum dots connected to ferromagnetic leads (at least 
qualitatively), one may use the decoupling scheme introduced in Ref. [135, 136 , but 
generalized by a self-consistent determination of the level energy to account for the 
exchange field 1 .'!!) 233 . Alternatively, one may calculate the exchange field first and 
then include it in the formalism by hand 234 . 

Figure[36]shows the differential conductance as a function of the transport voltage. 
For zero magnetic field there is a pronounced splitting of the peak in the parallel 
configuration [Fig. I36f a)]. which can be tuned away by an appropriate external 
magnetic field [Fig. RWb)]. In the antiparallel configuration, in turn, the opposite 
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Figure 36. Differential conductance in the P (a,d) and AP (b,e) configurations, 
and the corresponding TMR (c,f). Contributions to the conductance from spin- 
up (dashed lines) and spin-down (dotted lines) channels are also indicated. Left 
column (a,b,c) corresponds to zero magnetic field, whereas the right one (d,e,f) 
to the case when the compensating magnetic field B = 14.53 T is applied (for the 
Lande factor g = 0.152). The numerical results are for large U limit, while the 
other parameters are: k B T/F = 0.01, p = 0.2, T/D = 0.001. (After Ref. [233] ') 



happens, i.e. there is no splitting at B = [Fig. l3"6T d)] and a finite splitting at B > 
[Fig. I36f e)] with an additional asymmetry in the peak amplitudes as a function of the 
bias voltage. Figure RWc.f ) shows the corresponding TMR, which becomes negative 
without magnetic field and changes sign to positive when a compensating magnetic 
field is applied. 

Within the EOM approach the effect of spin-dependent quantum charge 
fluctuations is accounted in the self-consistent but intuitive manner. The real-time 
diagrammatic technique [1041 1167) enables one to construct a systematic approach, 
where the effect of ferromagnetic electrodes can be analyzed without any additional 
assumptions. Recently the resonant tunnelling approximation (RTA) was extended 
by Utsumi et al. [142 to account for influence of the electrodes' ferromagnetism on 
the Kondo phenomenon. This technique gives more reasonable results and allows for 
further systematic insight into the physics of the transport through quantum dots in 
the Kondo regime. 

The Kondo phenomenon in noncollincar magnetic configuration was considered 
in Refs |143l 11441 I145j . The results show that in symmetric systems, the Kondo 
anomaly gradually disappears when the magnetic configuration varies from antiparallel 
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to parallel one. However, the main drop of the Kondo peak appears already at small 
deviations from the antiparallel alignment. 

Recently Pasupathy et al. |5T| studied electronic transport through a single Cgo 
molecule attached to ferromagnetic Nickel electrodes in the Kondo regime. It was 
shown that the Kondo correlations appear even in the presence of ferromagnetic 
leads. The zero-bias anomaly in the noncquilibrium conductance, however, was split 
for the parallel alignment of the leads magnetization in agreement with theoretical 
predictions. For the antiparallel alignment, on the other hand, no splitting of the 
zero-bias anomaly was observed. Some residual splitting observed in this geometry for 
some samples can be interpreted as an effect of asymmetric coupling ^ T^. Similar 
behavior has also been observed in electronic transport through carbon nanotubes 
coupled to ferromagnetic leads |243j . 

7. Concluding remarks 

In this review we addressed certain aspects of electronic transport in mesoscopic tunnel 
junctions consisting of magnetic metallic nanoparticles, semiconducting quantum dots, 
or molecules attached to ferromagnetic electrodes. We are aware that the review does 
not cover all the problems considered in such systems and many references have not 
been included. There are some other aspects of electronic and spin transport which 
have been omitted. One of such aspects is spin torque acting on the island and/or 
electrodes in ferromagnetic single electron transistors based on magnetic metallic 
nanoparticles |244j . Another topic omitted in this review is the role of electron 
coupling to a phonon bath. Such a coupling leads to new interesting phenomena 
in transport characteristics [245]. For instance, the electron-phonon coupling leads to 
phonon satellite peaks in the differential conductance. However, the most important 
spin effects generated by spin-dependent tunnelling through the barriers separating 
the central part from the ferromagnetic electrodes have been addressed. 

We note that spin polarized transport through nanostructures, where the 
Coulomb effects become important, is rapidly developing in recent years. This 
applies not only to theoretical part of the subject, but also to the experimental side, 
where recent progress in nanotechnology allows to attach external leads to individual 
nanoparticles, quantum dots, and single molecules. An interesting and new field within 
this topic is the spintronics based on molecular magnets. Such molecules are good 
candidates for magnetic memory cells and also are considered as potential candidates 
for information processing devices. There is no doubt, that the rapid progress is the 
physics and technology of spin polarized electronic transport through nanostructures 
is stimulated by application possibilities in magnetoelectronics, spintronics, and 
information technology. 
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